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1.0 PROGRAM OBJECTIVES AND OVERVIEW 


The objective of the Propellant Nonlinear Constitutive Theory Extension 
Program is to develop and demonstrate an accurate and usable nonlinear thermal- 
mechanical constitutive law for solid rocket propellants. The program has 
been conducted through three phases of a four phase project. These four 
phases are summarized below: 


Phase I - Preliminary Study 


A variety of nonlinear theories were considered and five methods were 
selected for further study. 


Modified Swanson Theory 
Russian - (Hufferd) Theory 
Schapery Theory 

Gurtin Theory 

Quinlan Theory. 


Dr. Richard Wool presented a review of available micro-mechanics theories 


which could be considered for inclusion within the nonlinear constitutive 
theories. 


Phase II - Uniaxial Isothermal Investigation 

A series of unaxial tests were conducted with two materials - a PBAN and an 
HTPB propellant. The data was fit to each of the nonlinear theories. Then the 
material constants derived for each analytic method were used in a predictive 
calculation of a complex laboratory test history which included same typical 
long time rocket motor mechanical bonding sequences. A review of the analytic 
methods and predictive results showed that all of the theories could be adapted 
to give better correlations with realistic motor loading conditions. Other test 
histories were also suggested that would permit more direct evaluation of the 


pertinent nonlinear material behavior required for each theory. Numerical 
difficulties were encountered with same theories and refinements were determined 


that would hopefully eliminate these problems. All five theories were 
considered acceptable and were to be further developed in the next phase. 


Phase III - Two-Dimensional and Variable Temperature Investigation 

The analytic and numerical improvements were made with each theory. These 
refined methods were then evaluated with the test data from phase II plus the 
additional recommended tests from the phase II review meetings. Most of the 
analytic methods looked very good with the final refinements. This effort, 
using the latest data, used most of the scheduled time for this phase. In addi- 
tion to the suggested wmiaxial mechanical tests, a variety of shear, biaxial and 
uniaxial (some with simultaneous thermal and mechanical loading) tests were con- 
ducted with some very complex load histories. This latest test data was only 
considered in a preliminary way since the theory modifications were very 
exstensive. 


Phase IV - Three Dimensional Investigation (Not started) 

The additional laboratory data developed in phase III plus some 
instrumented structural test vehicle data will be used with the various 
nonlinear constitutive theories. One or more selected theories will be 
determined and the computer codes will be installed and checked out on the 
Air Force Rocket Propulsion Laboratory's computer. 


The program logic chart relating each of the phases and their respective 
tasks is presented in Figure 1. 
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y 2.0 TASK DESCRIPTION 

x: 2.1 PHASE I - PRELIMINARY STUDY 

= The objective of phase I was to critique the Quinlan theory, propose at 

¥% least one alternate approach to the constitutive law solution, and make detailed 
‘s research plans for evaluating and modifying the candidate constitutive law 

ae approaches. 

= 2.2 PHASE IT - UNIAXIAL/ISOTHERMAL INVESTIGATION 

ian The objective of phase II was to carry out modifications to Quinlan's law 
m, and to do theoretical development of the other candidate constitutive laws. 

x These laws were used to make stress-time predictions for uniaxial/isothermal 

i tests. Concurrently, actual uniaxial/isothermal tests were conducted in the 
laboratory. A comparison of the predictions and actual data was then made, and 
Hes the theoretical and experimental approaches were evaluated. 

a 

- 2.3 PHASE III - TWO-DIMENSIONAL AND VARIABLE TEMPERATURE INVESTIGATION 

The objective of phase III was to extend the theoretical development of the 
Pi candidate constitutive laws to two-dimensional and variable temperature tests. 
5 Stress-time predictions were made concurrently with actual testing in the labora- 
5 tory. A comparison of predictions and actual data was made and two constitutive 
law candidates were selected for further investigation from those results. 

: 

an 2.4 PHASE IV - THREE-DIMENSIONAL INVESTIGATION 

rs The objective of the final phase is to extend the theoretical development 
i of the candidate constitutive laws to three-dimensional and variable temperature 
a history tests. Concurrently, three-dimensional, variable temperature laboratory 
Bs tests are to be conducted. Data and predictions will be compared and any final 
os modifications to the constitutive law made. A validation experiment is to be 
conducted, and stress-time predictions made with the finalized version of the 

3 nonlinear constitutive law. A final assessment of the overall success of the 

= new law will be made and recommendations will be presented for further avenues 
= of research. Recommendations for utilizing the new law in existing solid rocket 
o motor structural analysis techniques will also be made. 
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3.0 LABORATORY TESTS, RESULTS, AND SUMMARY 


The laboratory testing was divided into the categories of uniaxial/ 
isothermal, two-dimensional and variable temperature, and three-dimensional 
investigations. 


The two propellants selected for the program were (1) a PBAN currently 
being used in the first stage of the Titan missile system (UTP~3001-750/7768) 
and (2) a HTPB propellant developed for the IUS motor (UTP-19,360B-400/1777). 
The first numbers are the propellant designation, the next the mixer size, and 
the last a batch number. 


In each of the uniaxial or biaxial groups, a specific test of each type was 
selected to show the test details in graphic and tabular form. While tests were 
run on both propellants, only one is shown. The details of test temperatures 
and rates are discussed with each test type. 


3.1 UNIAXIAL/ISOTHERMAL INVESTIGATION 

Testing wiaxial specimens of UTP-3001 and UTP-19,360B propellants, in 
phase II of the contract, was done for the nondamaged material as indicated in 
Figure 2 and for damaged material per Figures 3 and 4. Most of the tests were 
run with 1/2 x 1/2 x 6-in. bars with redwood end tabs. The exceptions were 
stress endurance (test 2) and constant rate (comparison to test No. 4) data were 
obtained with JANNAF Class B specimens. (1) The details of the individual test 
types are discussed in subsections below. 


The uniaxial bars were machined from redwood boxes of propellant. The 
redwood was sealed then lined in the same manner as a rocket motor. After a 
partial cure of the liner, propellant was cast in the box amd the system cured 
to provide a good bond to the redwood end tabs. The redwood box assembly and 
finished specimen are shown in Figure 5. After the specimen is mill finished, 


Reference 1 - Solid Propellant Mechanical Behavior Manual, CPIA Publication 
Wo. 21, Section 4.3.2 
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ws 24406R1 
r.. a 1/8 in. hole is drilled on the center line of each end tab for attachment 
ss purposes. (Attachment fixtures are shown in Figure 16). 
om Due to the enormous test load compared to the available time, both calendar 
rs and contract wise, the decision was made to run all tests on six-channel test- 
4 ers. Both the Chemical Systems Division (CSD) manufactured six-channel tester 
7, 
<< and a modified Instron were used in conjunction with a Hewlett-Packard computer 
sa: to collect digitized data for the tests. (See Appendix A details on the 
a automated data reduction system). While the CSD tester was equipped with an 
me oscillograph as backup for relatively short duration tests, the modified Instron 
a had no backup. In instances where a power surge occurred, the data were lost 
Na and the test had to be rerun. 
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Figure 5. Uniaxial Bar Specimen 
28804 

The attachment linkages on both testers were such that the specimens could 
not be put into compression when the crosshead was returned to an equivalent 
zero strain position. The strain measurement was done with a linear potentiome- 
ter attached to the crosshead; consequently, the data had to be modified to 
reflect the propellant strain relaxation behavior after the stress had returned 
to zero (free hanging specimen). 
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Strain relaxation was measured on samples in some of the tests during the 
final unload cycle. Cathetometer measurements were made periodically and strain 
versus time data were plotted. These data were used to estimate the relaxation 
behavior on cyclic tests where there was insufficient time for measurements. 


A data modification was made to estimate the peak or minimum stress and 
strain points, which were not recorded by the digitized data acquisition 
system. The sampling rate limited the crosshead rate that could be used 
and still obtain enough points to adequately define a ramp. The available 
computer memory also influenced the sampling rate in some of the long tests. 


3.1.1 Constant Rate Test No. 1 

Uniaxial constant rate tests to failure were rum on 6-in. bars of UTP-3001 
and UTP-19,360B. The 70°F tests were at crosshead rates of 10, 1, 0.1, 0.01, 
and 0.001 in./min. while 40 and 120°F tests were at 10, 1, and 0.1 in./min. A 
typical load-time curve is shown in Figure 6 for the UTP-3001 at 75°F and 
10 in./min., and tabular data are given in Table 1. Because of the computer 


200 


A- Stress (1/2 in. bar samples) Crosshead rate = 10 in./min. 
Test temperature — 75 F 
160 RAZ sera, Val m A 
120 


Stress. psi 


0 0.02 0.04 0.06 0.08 0.10 0.12 0.14 0.16 0.18 0.20 0.22 0.24 
Time, min. 


0 5.0 10.0 15.0 20.0 25.0 30.0 34.9 
Strain, % 


Figure 6. Test No. 1 = Straining to Failure for UTP-3001 750/7768 cs 
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data sampling rate, the 10-in./min. tests were run oe propellant at a time. 
Sample 4 (Table 1) shows the load cell reached the limit of its adjustment so 


{ did not record the specimen failure. The 6-in. bar specimens always fail below 
eS what would be obtained from JANNAF specimens. Since response properties are 

ee what is of interest, particularly in the small strain region, the uniform cross 
sectional area specimen does what it is supposed to do. The continual changing 


effective gage length of the JANNAF dogbone is avoided. 


These data were also reduced to secant modulus (\a/e) as shown in Figure 7. 


wt The data are compared to the stress relaxation modulus from test No. 4 later. 
oe 

vane 3.1.2 Uniaxial Stress Endurance Test No. 2 

un Stress endurance tests were run on the two propellants using JANNAF Class 
on B dogbones. The ultimate failure properties were of importance in this test 


rather than small strain response hence the dogbones. This is a constant load 
test with plastic extensometer to monitor the strain increase with time (also 
ve known as a creep test). They were run at 23, 70, and 120°F. The data shifted 


2,000 
A- Secant modulus (1/2 in. par samples) Crosshead rate = 10 in./min. 
a) Test temperature = 75 F 
oo 1,600 = 
As “ 
% F 
u wi 1,200 
2 
= 
: e 
: r Baas 
B 
400 
- 0 
- 0 0.02 0.04 0.06 0.08 010 0.12 O14 016 018 020 0.22 0.24 
7 Time, min. 
a, 
—_ 0 5.0 10.0 15.0 20.0 25.0 30.0 34.9 
oe Strain. % 
- Figure 7. Test No. 1 ~ Secant Modulus for UTP-3001 750/7768 
ss 28751 
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to a master endurance curve as shown in Figure 8 for UTP-19,360B with typical 
70°F data given in Table 2. The strain creep effect on stress is shown in 
Figure 9 where the engineering stress (F/A,) is multiplied by the extension 
ratio (1+€) to account for the sample's necking down. The secant modulus (Ag/€) 
for the 70°F data are shown in Figure 10 and data are given in Table 3. This 
same type of data were generated for UTP-3001. 


orre49r 
nn CA 


3.1.3 Multirate Test No. 3 

Constant rate tests were run on the two propellants and the rate was 
changed in the middle of the test. Because of the different response for the 
high-low compared to the low-high, an example of both is given. The test was 
included with the nondamaged tests because there was no rest or reversal in the 
crosshead direction. The first leg of the test could be considered the damage. 
The 1.0 to 0.1 in./min. rate change is shown in Figure 11 for UTP-19,360B and 
the 0.1 to 1.0 in./min. is shown in Figure 12. The corresponding data for the 
first sample of each group are given in Table 4 and 5, respectively. 


3.1.4 Stress Relaxation Modulus Test No. 4 

The stress relaxation modulus tests were run at a nominal 3% strain using 
1/2 x 1/2 x 6-in. samples of propellant bonded to redwood end tabs for both 
propellants. The samples were loaded to 3% strain at a crosshead rate of 
1 in./min. for temperatures of 20, 43, 73, and 122°F. The load was monitored 
with time. Strain was determined by cathetometer measurements on the samples. 
The relaxation modulus (Ep) was then calculated at specific time intervals by 
the following equation: 


F 
ER =— A 
Ag € 


where F = force as measured by the load cell 


A, = initial area 
1+e 


strain 


» 
" 


The master stress relaxation modulus data for UTP-3001 are presented in 
Figure 13 and typical data at 73°F are given in Table 6. 
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iy TABLE 3. TEST NO. 2 - UNIAXIAL STRESS ENDURANCE CREEP BEHAVIOR FOR 
< UTP-19, 360B-400/1777 at 70°F 
ca 


(SHEET 1 OF 2) 


T8706 
Secant 
Temper- Stress Modulus — 

oe ature, Time, Strain, Load, Area, (F/A), do, (do/e), 
F Sample pin. $ g in.@ psi psi psi 
of 70 1 0.13 20 7,180 0.185 85.5 102.6 513 
“7 0.33 30 111.2 370 
re 0.87 40 119.7 299 
2.0 50 128.3 257 
ay 60 136.8 228 
By 5.0 61 137.7 226 
: 70 2 0.27 50 8,720 0.184 104.4 156 .6 313 
eos 0.41 60 167.1 279 
ais 0.55 64 171.2 268 
La 70 3 1.0 30 6,000 0.179 73.8 95.9 320 
: 5.3 40 103.3 258 
= 11.0 46 107.7 234 
. 14.9 50 110.7 = 221 
“ 21.0 50 110.7 221 
% 70 4 7.5 33 5,200 0.182 62.9 83.6 253 
ae 18.0 38 86.8 228 
! 95.0 kg 93.7 191 
~ 196.0 56 98.1 175 
xm 247 .0 56 98.1 175 
Ped 70 5 65.0 27 4,480 0.180 54.8 69.6 258 
4 190.0 30 71.2 237 
f 440.0 35 73.9 211 
ae 1,640.0 40 76.7 192 
SS 3,325.0 45 79.5 177 
6,086.0 47 80.6 172 
oy 70 6 2.2 21 4,750 0.177 59.1 71.5 340 
4) 14.0 29 76.2 263 
Ns 30.0 32 78.0 244 
ry 84.0 37 80.9 219 
_ 150.0 38 81.6 214 
: 833.0 50 88.7 177 
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MS 
ae TABLE 3. TEST NO. 2 = UNIAXIAL STRESS ENDURANCE CREEP BEHAVIOR FOR 
= UTP-19, 360B-400/1777 at 70°F 
4 (SHEET 2 OF 2) 
= 78706 
25 
.%: 
oN Secant 
fa Temper- Stress Modulus 

ature, Time, Strain, Load, Area, (F/A), dAo, (AG/E), 
. F Sample ‘min. % g in.2 psi psi psi 
= 70 7 0.7 25 5,630 0.183 67.8 84.8 339 
= 4.0 30 88.1 294 
=i 10.0 35 91.5 261 
‘i 26.0 42 96.3 229 
b.. 92.0 51 102.4 200 
ae 177.0 57 106.5 187 
ee 70 8 10.0 24 4,660 0.184 55.8 69.2 288 
os 77.0 30 72.5 2he 

140.0 32 73.7 230 

io 1,418.0 yy 80.4 182 
a: 2,900.0 49 83.1 169 
= 3,500.0 49 83.1 169 

70 9 10.0 19 4,030 0.178 49.9 59.4 312 

30.0 20 59.9 299 
at 70.0 22 60.9 277 
Po 5,496.0 35 67.4 193 
S 10,040 37 68.4 185 
3 Unbroken 
28 ,600 me 4 70.4 172 

ro 3.1.4.1 Constant Rate Modulus 
" Constant rate modulus tests were run on the 6-in. bar samples described 
2 above and on JANNAF Class B specimens. The 6-in. bar samples use the wood-to- 
“ wood distance as the effective gage length while the JANNAF's were analyzed 
a using 2.70-in. effective gage length through the test even though it is varying. 
> Strain was determined by the crosshead travel. The constant rate modulus (F(t)) 
“et is calculated from the following equation at specific strain levels through the 
tests. 
r,s F r 
Po! F(t) =<— — 
x. Ay €(t) 
oj where: ¢(t) = strain at time t 
om 
“, 
e 
ee 20 
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Test No. 3. High- 
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UTP.19, 360B~400/ 


Low Co ate Tests of 
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1.5 Multiple Loading Test No, 5 


The Multiple loading tests on 6-in, 


Period between cycles, 
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this data had to be 


@ 
(-% 
2 
~ 
o 
— 
n 
7) 
© 
= 
= 
wn” 


Strain, % 


12. Test No. 3 - Low-High Constant Rate Tests of 
UTP-19, 360B~400/1777 at 70°F 


28392 


The rework of the data consisted of estimating maximum and minimum stress 
and strain points that were not picked up by the digitized data acquisition 
system. The data reduction system computed strain from crosshead travel. This 
was satisfactory except at and below zero stress. The sample linkage attachment 
was such that specimens would not be put into compression with the exception of 
test No. 11, which was run in an Instron with rigid clamp jaws. The actual pro- 
pellant strain decay was estimated from other tests where strain recovery was 
monitored by cathetometer measurements for the part of the tests at zero stress 
(1.e., no load on the samples). The 5-in./min. crosshead test for UTP-19,360B 
was selected as typical and shown in Figure 15. While the load-unload ramps 


TABLE 4. TEST NO. 3 - HIGH-LOW CONSTANT RATE TESTS OF UTP-19,360B—400/1777 
AT 70°F WITH 1/2 x 1/2 x 6-IN. BAR SAMPLES 
T8707 


Sample No. 1 
Crosshead speed, in./min. 1.0 and 0.1 
Chart speed, in./min. 10 and 1 


Load scale 5 lb/in. 
Area, in. 0.500 x 0.501 
Gage length, in. 


Chart Distance, Stress (F/A,), 
in. % psi 


0 
5.190 
10.379 
19.960 
23.952 
27.745 
31.337 
34.731 
37.924 
40.719 
4.15 2.10 41.916 
Change to 0.1 in./min. 

1.74 34.731 
1.70 33.932 
1.68 33.533 
1.66 33.134 
1.65 32.934 
1.65 32.934 
1.73 34.531 
1.83 36.527 
1.93 38.523 
2.12 11.667 42.315 
2.29 13.333 45.709 
2.43 15.00 48.503 
2.55 16.667 50.898 
2.66 18.33 53.094 
2.74 20.00 54.691 
2.76 70.50 55.090 
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TABLE 5. TEST NO. 3 - LOW-HIGH CONSTANT RATE TESTS OF UTP-19,360B-400/1777 


AT 70°F WITH 1/2 x 1/2 x 6-IN. BAR SAMPLES 


Sample No. 4 
Crosshead speed, in./min. 
Chart speed, in./min. 


Load scale 5 1b/in. 
Area, in. 0.502 x 0.500 
Gage length 6.00 


Chart Distance, Load Signal, Strain, Stress (F/A)), 
in. in. t psi 


0 
2.789 
5.578 
8.566 
11.355 
13.745 
16.534 
18.725 
21.315 
23.705 
25 896 
27.092 


33.466 
40.837 
4H U22 
48.805 
51.195 
56.175 
60.359 
15.634 63.347 
17.300 65.737 
18.967 67.530 
20.634 69.323 
22.300 70.717 
23.967 72.112 
25 634 73.307 
27.80 74.701 


r>erreseeseeosso 


e 
WW OWOWOAUEN — 
Hoon Fworw we 


@ 


e 


nge to 1 in./min. 
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0 
0.4 
0.8 
1.2 
1.6 
2.0 
2.4 
2.8 
3.2 
3.6 
4.0 
4.1 
Cha 
4.2 
4.25 
4.3 
4.3 
Ky 
4.5 
4.6 
4.7 
4.8 
4.9 
5.0 
5.1 
5.2 
5.3 
5.4 
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10.000 
Nominai strain = 3 0% 


6.000 Temperature F Log A, 


Oo 20 
© 43 
A 73 
© 122 


3.000 


1.000 


Stress relaxation modulus 
(Eg). psi 


V/A, min 


Figure 13. Master Modulus Data for UTP-3001-750/7768 with 
Experimental Shift 
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are not clearly shown in the figure because of the scale necessary to show 
the total test duration, the detailed data are given in Table 8. 


3.1.6 Creep Test No. 6 

Creep tests were run on 6-in. bars of UTP-3001 and UTP-19,360B by hanging 
weights on the samples and allowing them to deform under load. The loading was 
done by setting the weight on the Instron crosshead and then lowering the cross- 
head away from them as shown in Figure 16. Proper spacing of the segmented 
weight provided a means of incrementally unloading the samples. The two propel- 
lants were run separately to allow access for spacer insertion without disturb- 
ing the samples. The procedure for unloading part of the load and still 
maintaining the balance as a constant load creep test is shown in Figures 17. 
Tests were run at the maximum load the sample would be sure to survive, then at 
half and quarter loads. Tests were repeated for 70, 120, and 40° F. The full 
stress creep test for UTP-19,360B at 71°F is shown as typical. The strain time 
curve from crosshead and cathetometer measurements is shown in Figure 18 with 
engineering stress, corrected stress, and secant modulus shown in Figures 19, 
20, and 21, respectively. The test data are given in Table 9 and summarized in 
Table 10 with cathetometer and secant modulus data. 

(Text continued on page 29.) 
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MASTER MODULUS CURVES FOR UTP-3001 STRESS RELAXATION AND 


CONSTANT RATE TESTS 
(SHEET 1 OF 2) 


6-in. Bar, 
ER, psi 


8,750# 
7,600* 
6,850# 


6,350* 
4 ,000# 
3,500 
3,170 


3,000 
2,400 
1,930 
1,720 
1,570 


1,470 
1,220 
1,020 
910 
840 


800 
675 
580 
520 
490 


470 
412 
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6~-in. Bar, 
F(t), psi 


T8709 


JANNAF , 
F(t), psi 
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De TABLE 7. MASTER MODULUS CURVES FOR UTP-3001 STRESS RELAXATION AND 
ay CONSTANT RATE TESTS 


Ce (SHEET 2 OF 2) 
ro 18709 | 
om Reduced Time, 6-in. Bar, 6-in. Bar, JANNAF , : 
a t/Ap, min. Er, psi F(t), psi F(t), psi 
: 1x 102 230 353 332 
: 2x 102 218 326 310 
a 4 x 102 207 300 290 
5 6 x 102 201 290 280# 
8 x 10° 198 280 273* ) 
“a 1x 103 193 275 268" | 
c 2x 103 185 2588 2558 | 
eS 4 x 103 182 245@ 242% | 
aS 6 x 103 179 236# 235* 
we 8 x 103 172 2328 2308 : 
1x 104 169 2288 2278 | 
Ms 2x 104 163 2198 2188 | 
7 4 x 104 1578 2108 2108 | 
x 6 x 104 1548 205# 205* 
sh 8 x 104 1528 203# 202# ! 
r 1 x 10° 151# 2008 200* : 
a 2 x 10° 1478 1938 1948 | 
a 4 x 109 1438 1888 1888 | 
be 6 x 10° 1414 184# 185* 
e 8 x 10° 140# 1828 183# 
“ 1x 106 139# 181# 182# 
: 2x 10° 1368 1778 1738 
RN: 4 x 106 133* 1748 1748 
a 6 x 106 1318 1718 1728 | 
8 x 106 130 169® 170# | 
+, ® Extrapolated data 
A 3.1.7 Cyclic Loading Test No. 7 


a 


Cyclic loading tests were run on 6-in. bars of UTP-3001 and UTP-19,360B pro- 
pellants at ambient temperature. The cycling was for 20 cycles at nominal 


Fos ’ 


iia 
. 


strain levels of 4, 8, and 12% for UTP-19,360B with UTP~-3001 limited to 12%. At 
pt the end of the test (after unloading to zero stress) the strain was monitored on 
a the samples with a cathetometer. The test at a nominal 12% strain is shown in 
Figure 22 for UTP-19,360B. These data were modified to insert the estimated 

“ maximum stress points and the propellant strain decay while at zero stress. 

a Data for the test are given in Table 11. (Text continued on page 34.) 
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TEST NO. 5 - 1/2-IN. BAR STRESS WHILE CYCLING 
(SHEET 2 of 3) 
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TEST NO. 5 - 1/2-IN. BAR STRESS WHILE CYCLING 


TABLE 8. 


(SHEET 3 of 3) 


SSS COGS OVO € NOOO ict OF wt {IO O28 OS OS 28 OSES SHOP ROODNOMSOIIS NS OS SSSOOHSEC SoS 
SOS SOsSoNank JMNM BNA OT SOOO SOSSS SOSOSC SH UWAMOROVMMADATOOo cooecooesco 
O00 OO COOK CE ROT EM CWA O8 OO 09 O90 O08 OD EE NE OULU Rw HOSS OOOO OO80S60 


Cs i rr? 


O80 OS S68 OO BS OS ewi ef NOSSO SO &O OS OS SS SCO SSe Soe oo Sconanoscoss Sooecocoeooe 


SES OCSCS OOM AUILAN NR CONICS OS 96 SOS 8OSOS OF OH DANA NO OM AMO SOCOCSCOCeCSe Sco 

S66 ©OSSS OA CSMMIUG OV ON CO WO O80 20 88 OS OOSS CANN MEM AMNMAR RS OOOSC OS Oc OO 

ee ee a OO en a Be ee 
MUNMTNOOEM NM EMNOR Meet 


SSS CoCeOCMNeN Dooce to [fon SoCo Sole SOCeOSCEIMMNNOOT OMe Dow Ooocecscescesecs 
O90 E90 OE NL DS FO. 09 NA ND UD OS 90 00 06 0896 CNNOD RO ODODE AV. OOOS OS 5000 


SOS COSCON A OM HON NO DO IM SOSS SSCS ESOS OPIN DE MIDONO Ne Ssessescos 
MOICMEN OS Cit MINIM TRON Mate 


S02 SSSCOCNSCamMVOTMSSUNNOUNSCSSCSSSSSe OSS Oma NNNYOSeCVeMess ocesooosscs 
BOS OO9S SVEN WOMAN END ARS OOSOSOOS COOOONNN DOMDN we ODESS COCs oSoScS 


Seo ScoC Onda OMe HCO ROSIN OS SS OSes COSC ONINNOUINS LOO nmnme cooSoSScco] 
AN NMOS Mews QU CMEL? CHM tet 


pentapeptide echendy tech apebeptemrdoe ns gab nb ded od ly opamp 
Soc cooeeof On ¢€-Oy.KW ONT DINO OS OO 9SSS 9886 CNOV CNOS DOMAIN AS OBOOCSSSSS 


SOSSSSOOMMNSY MOSK O NOUN SSSSSOSSSSOSONE HRMSN TOO MING SSSSSsSseS 
CAN DSP LS LP: PANURSMICMN GNM tet 


SS CUNO MR ORSS C OS CUO. CLASS EMME EUAN IDS SODANDS OURS. FBO. Oe OU 
MEM OSOS CUNG NENA MOM MAD OM AN MAMMA SOSSS TOM AOMRMR OC OM Nt ONOM OMMw OS 


a ee Cr ee 


SSSSOSSOSA FLINOE POD OUEMNNNSSSSSSSSS SMM LOO ANNSONE MMinasocse 


BOK MONE WNL NSS BOOTHE 
PIN TMM STS weeeeee a in inal ced Sat Sal 
ARAARARL AARARRAR ARSARABRAN 
mON STING FOCOON OH sonoonm 

[nt a ad MMM Mmmm mentee 
RRRRRRR RAARARARKR AAABRAARK 


pS et Rated ty Sut ay $2 Re te pA pat de be Peat AT Ne Rie fleck cohol ap lla Ray ie Bef ely to ed 
Mat DMMSMONRD NERONAMNONAMNO AME NUNDUORSUAMe OmORODONR ON Di OVEROMAUON OF 
ROMMe hoM@ma moMoM BO ei Sn OR OO em de eh Dont Jie Ono 
BM tate IRDSNMMS NDS amen ORR SARROCOHMUR DONE OR DON TO DOMANNANMMN 
SPONGES MAN MACS SO00ORR BRARS PUNO MDMNMMMMMSe eee eT OOO TNS eo 

ORONO ATLA ALLAN LAUT AUT LALLY OR UD AAU LEAL LALLA LLALALA ON 00,60 
BSGGRADASADDS ROPOARA KOK PARRA AAAOAOO OONINNANN RUN OU CUCU CUCU RUCUCUCG RUCUEURUCUCUEUS Oe 
AAAAnNAMA AoA Melee Citric etoiwti net 


NOMCNOKRDOOANMENORE COANMENORDH OAM ENOR DOS AVM EMNORDHORUMEINOKR DKON 
KRARRARRAOWD DDDOOS DOO. 20.0 0-0 OO BOOS OOOCS OBOE at ot whet ot or wt ot HIN CUNION UCI AES PPS 
Add Ade Hae el tet el eletet tate el etetet ei et 


33 


3.1.8 Relaxation Test No. 8 

Stress relaxation tests were run 
on 6-in. bars of UTP-3001 and UTP- 
19,360B propellant for a 24-hr period 
and then monitored for strain decay 
after they were unloaded. The tests at 
ambient temperature were repeated for 
4, 8, and 12% nominal strain levels for 
UTP-19,360B but limited to 4% for UTP- 
3001. They were loaded at a crosshead 
rate of 20 in./min. and wloaded at 
1 in./min. after the 24-hr relaxation. 
A plot for the UTP-19,360B at 4% nominal 
strain is shown as typical in Fig- 
ure 23, and data are given in Table 12. 


3.1.9 Predamaged-Relaxation Test No. 9 
The predamage-relaxation tests 
were run with 6-in. bars on UTP-19,360B 

propellant at ambient temperature. 
Figure 16. Test No. 6 - Creep They were preloaded to 12% and unloaded 
Pony wath: O30) Par sce taaaerste at a crosshead rate of 0.1 in./min., 

allowed to rest, then reloaded to 8 or 4% strain at 20 in./min. After relaxing 


1 hr the samples were wloaded at 1.0 in./min. and strain was monitored 


19065-1 


after unloading. These tests were repeated for a 6% predamage strain followed 
by a reload to 4 or 2% as above. The 12% predamage followed by 8% stress 
relaxation for UTP-19,360B is shown in Figure 24, and tabular data are given 
Zt in Table 13. Comparison data for UTP-3001 was limited to 6% predamage followed 
i by 3% relaxation. These data were modified to obtain the peak stress and 
strain relaxation after the samples were unloaded. Strain was monitored 
with a cathetometer after the relaxation part of the test. 


‘ 3.1.10 Complex Multiple Load Test No. 10 
The complex multiple load tests were run with 6-in. bars on UTP-3001 and 
UTP-19,360B propellants. Tests were run at crosshead rates of 5, 1, and 
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Figure 17. Procedure for Loading and Incrementally Unloading Creep Samples 
28805 

0.1 in./min. The test sequence was 12 to 8 to 12 to 4% strain, then unload, 

reload to 4% strain, and wload (four cycles) with cathetometer monitoring of 

strain decay. The same type of sequence was repeated with maximum strains of 

8 and then 4% where the 4% maximum strain was shortened by one cycle. The 

5 in./min., 12% maximum strain test on UTP-19,360B as typical is shown in Fig- 


ure 25, and tabular data are given in Table 14. The data were reworked to insert 


maximum and minimum stress values as well as strain decay for unloaded specimens. 


The cathetometer strain after the final unload was incorporated into the data. 
The 5 in./min 12% strain test for UTP-3001 was deleted because the modulus 
was outside the range of the remainder of the data (i.e. carton to carton 
difference). 
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Figure 20. Corrected Stress-Time Data for Creep Test of 
UTP-19, 360B-400/1777 at 71 F 
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3.1.11 Quinlan Complex History Test No. 11A 

This particular test is very complex and required 2 full months to com- 
plete. It was run with single samples (one at a time) in an Instron using 
6-in. bar specimens of UTP-3001 and UTP-19,360B propellant. The samples 
were pinned through the redwood end tabs and clamped onto pins to avoid crush- 
ing the redwood (because that would generate a compression load on the sample) 
during attachment. All coupling joints were heavy and tightly pinned so 
that the sample would be put into compression when returned to the zero strain 
position. 
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‘ Figure 21. Secant Modulus Data for Creep Test of 

UTP-19 , 360B-400/1777 at 71 

a 28752 
#8 The test on UTP-3001 is presented but the test is so complex that it 
pe. has been divided into three parts. Part 1 is described in Table 15 with 

; data given in Table 16. Since an actual Instron trace of the load-time curve 
2 was obtained, the peaks and minimum (compression) stresses were selected 

a data reduction points. Part 1 of the test is expanded in time scale to show 
re some detail of the process (Figures 26 through 29). Part 2 is described 

in Table 17 with data given in Table 18. Test sequence for Part 2 is shown 

na in Figures 30, 31, and 32. Part 3 (selected cycle maximum and minimum) are 
me given in Table 19. The last figure of this part (Figure 32) is some of the 
M eyclic loading at the end of the test. The chart speed was set such that 
9 good definition of the cycles could be recorded. 
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TEST NO. 6 - 1/2-IN. BAR STRESS WHILE CYCLING (FULL STRESS CREEP) 
(SHEET 6 OF 8) 
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a 
MS 
ce TABLE 15. TEST NO. 11, PART 1 - QUINLAN COMPLEX HISTORY FOR UTP-3001 
Rate, T7859 
hy: Cycle in./min. Remarks 
“st 7 enna 
SS 1-7 Load 2 Approximately 15 min. rest after cycle 
cee Unload 2 Approximately 15 min. rest after cycle 
zt 8 Load 1 Approximately 15 min. rest after cycle 
bo Unload 1 
5m 
ae 9 Load 5 Approximately 15 min. rest after cycle 
< Unload 5 
10 Load 0.5 Approximately 15 min. rest after cycle 
>s Unload 0.5 
ae 
nN 11 Load 10 Approximately 15 min. rest after cycle 
aN Unload 10 
om 
fe 12 Load 0.2 Approximately 15 min. rest after cycle 
ae Unload 0.2 
a 13, Load 2 
bt Relax 1/2 hr 
Unload 2 Approximately 30 min. rest after cycle 
5) 14 Load 2 
A Relax 1 hr 
a) Unload 2 4 day rest after cycle 
Ca) 
15 Load 2 
os Relax 1 hr 
i Unload 2 7 day rest after cycle 
pas 
"3 
pay The later part of the cycling is represented only by the maximum and minimum 
: stress-strain points. Part 3 of this test, the balance of cycling to failure, 
a 4s recorded in Table 19 as maxima and minima for selected cycles sufficiently 
Se close to describe the upper and lower bounds. A plot of the data would be 
similar to Figure 32. The strain values in Table 19 are stable, while the maxi- 


mum stress shows a continual decay, and compressive (negative values) stresses 
te, become less compressive. 
(Text continued on page 76) 
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TABLE 17. TEST NO. 11, PART 2 = QUINLAN COMPLEX HISTORY FOR pre pene 


r 714 
me Cycle Rate, in./min. Remarks 

"3 

- 16 Load 0.02 

% un load 0.02 Approximately 30 min. rest after cycle 

17 Load 0.02 
Es unload 0.02 Approximately 30 min. rest after cycle 

ss 18 Load 0.05 

“f relax 3 hr 
= unload 0.05 2 weeks rest after cycle 

= 19 Load 0.02 

ce unload 0.02 Approximately 30 min. rest after cycle 

> 

at 20 Load 0.02 

il unload 0.02 Approximately 30 min. rest after cycle 
se 21 Load 0.05 

at relax 3 hr 

: unload 0.05 1 month rest after cycle 
a 22-42 = Cycling 5 Several cycles monitored followed 
7 by several with only maximum and 

sy minimum recorded. 

%) 

fo 

> Note: Part 3 of this test was continuing cycling to failure. The maximum and 
¥ minimum values at larger intervals have been tabulated but not incorporated into 
* the disk files. 
x 3.1.12 Similitude Test No. 12 
a This similitude test and the one following it were run with the intent that 
rs] only the strain-time history would be supplied to the subcontractors, who would 


then predict the stress-time histories from their respective predictive theories 
for the program review meeting held at the end of the phase II testing in 


Bf? 
¥ 
.. 


pia) 

ix September 1982. 

is 

© This similitude test on the 6-in. bar specimens of UTP-3001 and UTP-19,360B 
N propellants was: (1) a 0.01 in./min. ramp to 10% strain followed by relaxation; 


‘ (2) a 1 in./min. unload to 5% strain and relaxation; and (3) 0.1 in./min. ramp 
to failure. The same test was repeated for ramp rates of 0.001, 0.1, and 


a 0.01 in./min. The data as reported at the September meeting are shown for UTP- 
i 19,360B in Figure 33. The ramp rates were 0.01, 1, and 0.1 in./min. Tabular 

s data for the test are given in Table 20. 
“N 
ty 
rr 
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TABLE 19. TEST NO. 11, PART 3 - UTP-3001-750/7768 1/2-IN. 
BAR STRESS WHILE CYCLING 


(SHEET 1 OF 2) 
T8715 
Time Strain Stress Remarks 

74388 .0479 2.40 -6.37 End of cycle 20 
74390 .3479 12.17 93.23 Peak of cycle 30 
74390 .4539 2.40 -5.98 End of cycle 30 
74391.3819 12.17 92.03 Peak of cycle 34 
743914239 2.40 -5 58 End of cycle 34 
74391.9432 12.17 91.24 Peak of cycle 40 
743919859 2.40 -5 .58 End of cycle 40 
74392 .0799 12.17 89.64 Peak of cycle 50 
74392 .1219 2.40 -5.18 End of cycle 50 
74393 .0099 12.17 88 45 Peak of cycle 60 
74393 .0519 2.40 -5.18 End of cycle 60 
74394 .7799 12.17 84 86 Peak of cycle 80 
74394 .8219 2.40 4.78 End of cycle 80 
74395 .7279 12.17 80.10 Peak of cycle 102 
74395 .7479 2.40 -4 78 End of cycle 102 
74418 63479 12.17 76.49 Peak of cycle 164 
74418 £3979 2.40 -4 .78 End of cycle 164 
74451 .5979 12.17 71.71 Peak of cycle 330 
744518779 2.40 -4 38 End of cycle 330 
74474 £2779 12.17 68.92 Peak of cycle 442 
74474 .3179 2.40 -3.98 End of cycle 442 
74698 .0179 12.17 62.95 Peak of cycle 980 

2.40 3.59 End of cycle 980 
74928 .0179 12.17 60.16 Peak of cycle 1980 

2.40 -3.59 End of cycle 1980 
75158 .0179 12.17 56.97 Peak of cycle 2980 
75388 .0179 12.17 56.18 Peak of cycle 3980 

2.40 -3.59 End of cycle 3980 
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2) 
‘] 
ne TABLE 19. TEST NO. 11, PART 3 - UTP-3001-750/7768 1/2-IN. 
‘ BAR STRESS WHILE CYCLING 
: (SHEET 2 OF 2) 
18715 
Time Strain Stress Remarks 
Pet 
75618 .0179 12.17 53.39 Peak of cycle 4980 
s 2.40 3.59 End of cycle 4980 
75695 .5179 12.17 51.39 Peak of cycle 5317 
% 2.40 -3.59 End of cycle 5317 
3 SAMPLE BROKE 
3.1.13 Three-Step Relaxation Test No. 13 
ce The three-step relaxation tests were run as similitude tests in which the 
strain-time data were reported to the subcontractors and they were to predict 
"4 
“4 the stress-time histories. The test consisted of loading 6-in. bar specimens at 


0.05 in./min. crosshead rate to 10%, relaxing 1 hr, unloading at the same rate 
‘ to 7% strain, relaxing 1 hr, and repeating the process at 3% strain, then 

3 unload. The tests were then repeated with nominal 24 hr relaxation periods for 
’ both UTP-3001 and UTP-19,360B propellants. 


These data were reworked to make sure that peak and minimum stress points 
: were included in the data. Relaxation was monitored after the sample was 

, unloaded to zero stress. Data for UTP-3001 is shown in Figure 34 with the 1 hr 
relaxation periods. Tabular data are given in Table 21. 


(Text continued on page 103.) 
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3.1.14 Propellant Aging Effects During Phase II Testing 

The numerous and complex tests involved in the uniaxial-isothermal 
evaluation of the two propellants (UTP-3001 and UTP-19,360B) covered an approxi- 
mate 3-month time period. Since the 6-in. bar specimens at constant rate (test 


No. 1) were run first, they were used as the standard to compare the balance of 
the data. 


Stress-strain plots of UTP-3001, test No. 1 are shown in Figures 35, 36, 
and 37 for 41, 75, and 124°F tests, respectively. Crosshead rates of 10, 1, and 
0.1 in./min. are shown for each temperature with the ambient tests down to 0.001 
in./min. The same data are shown in Figures 38, 39, and 40 for UTP-19,360B. 
Comparison of the other tests was based on the initial ramp loading or undamaged 
state. Test No. 6 was a full, one-half and one-quarter load creep test with 
crosshead loading rate of 1 in./min. The peak stress-strain points, when the 
crosshead stopped, were compared to the constant rate plots as shown in Figure 
35, etc. All the data are given in Tables 22 and 23 for UTP-3001 and UTP- 
19,360B, respectively. When crosshead rates did not match, a time-temperature 
equivalent stress value was selected for comparison purposes. These data were 
considered approximate. Stress values that were considered reliable (i.e., very 
close) were marked in the tables with a single approximate sign (~) while those 
that left some doubt due to the shift were marked with a double approximate (=). 
Some tests were run at 5 in./min. rather than the 10 in./min. comparison data. 
The 10 in./min. equivalent of the 5 in./min. data can be determined from the 
rate shift (i.e., log Ap = log 10/5 = 0.30). The WLF curve (in text figure) is 
used to pick off the equivalent temperature. The high rate stress value being 
sought is equivalent to a lower temperature at the lower rate (5 in./min.). 
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The equivalent temperature is then used to obtain the equivalent stress 
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value from a stress-temperature plot (in text Figure). 
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The stress-temperature plot can also be used to obtain equivalent stress 
values when the temperatures do not match at the same rate. 
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After the initial ramp loading, the propellant was considered damaged. In 
most instances the tests were of short enough duration to neglect any aging 
effect during test. Since bulk storage at controlled ambient conditions has not 
shown a significant aging effect over 3 months, this was also neglected. Rather 
it was lumped into the between carton difference. The sample handling procedure 
established was to machine one carton at a time, hold all samples overnight in 
a nitrogen flushed dry box (<10% RH) before testing, and then test all samples 
from the dry box before machining additional specimens. The data in Table 22 
for UTP-3001 indicate little difference between boxes 1 and 2 but a one-third 
higher average stress in box 3. The other significant change was within a box 
due to the amount of storage time in the dry box. During the 13 days the resid- 
ual of box 2 was in the dry box, the stress changed from 30% below the mean to 
10% above the mean. The UTP-3001 tests from box 3 have been replaced with 
another carton which was sample tested before doing the program tests. This 
provided a reasonable assurance that the results would be of the same family as 
the original constant rate data (test No. 1). The data from UTP-19,360B in 
Table 23 do not show as much effect for dry box storage or between box 
differences compared to the UTP-3001 propellant. 


Both propellants have indicated that dry box storage will increase the 
stress capability. This is assumed to be due to loss of moisture and it is 
reversible as is shown in Figure 41 for UTP-3001. In this example, specimens 
held in the dry box for 4 months were exposed to 58 and 85% relative humidity 
for 1 week. All six samples were tested simultaneously in a CSD multistation 
tester. 


The one test which may need some special treatment is the Quinlan complex 
history No. 11 (see section 3.1.11). After 12 cycles and 2 relaxation periods, 
the samples were removed for 4 days storage in a dry box. They were then 
followed by (1) relaxation - 7 days in a dry box, (2) two cycles ~ a relaxation 
period - 2 weeks in a dry box, (3) two cycles - a relaxation period - 1 month in 
a dry box, and (4) cycling to failure. At this point it is not clear just how 
to separate the stress increase due to rehealing from the dehydration effect of 
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Figure 35. Initial Ramp for UTP-3001-750/7768 39°F Tests Compared to the 
6-in. Bar Constant Rate Data 
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the dry box particularly for the 1 month storage period. The test times at 
approximately 50% RH were too short to have any significant effect on propellant 
stress capability. Both propellants showed reasonable correspondence of the 
initial ramp to the reference constant rate data. 


3.2 TWO-DIMENSIONAL AND VARIABLE TEMPERATURE INVESTIGATION 

The biaxial and nonisothermal testing was conducted on specimens of UTP- 
3001 and UTP-19,360B propellants as detailed in Figure 42. The biaxial samples 
were cast into prelined redwood boxes with a 1.25-in. gage length by 6-in. wide 
and machined flat to a 0.25-in. thickness. The response properties rather than 
failure properties were of interest so the discontinuity at the redwood inter- 
face did not affect the desired behavior. The 1/2 x 1/2 x 6-in. specimens were 
used for straining-cooling and cyclic strain tests. Shear relaxation tests 
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ra were run with 1 x 1 x 3-in. specimens bonded directly to steel anvils. Details 
oy are given in later sections. The equipment discussed in section 3.1 was 
aa utilized for this testing; however, only three biaxial specimens could be tested 
bs at once because of space limitations in the oven. 
ea The biaxial specimens used in this part of the program were cast into 
a redwood boxes similar to that shown in Figure 5. The space between redwood 
f blocks was 1.25-in. instead of the 6-in. for the uniaxial bars. A mill finished 
a specimen is shown in Figure 43. The propellant was left flat rather than 
ne necking it down as with standard JANNAF biaxial specimens. The gage length was 
Py designated at the wood to wood distance for strain evaluation. Response 
i 


properties rather than failure properties were of interest so the boundary 
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ef 

perturbation was neglected. Some data taken from the literature (2) have been 
. reproduced for evaluation of the stress-strain behavior across the biaxial 
g field (see Figures 44 through 47). 

The shear samples (Test No. 17) were 1x 1 x 3-in. blocks of propellant 

Po that were bonded to the test fixture (shown in Figure 48) after being 
a machined. The pull rods were attached to the offset plates so that the load 
A. 
—~ was transmitted through the center of the sample as shown. Since strain 
a) was limited to 5% for the shear relaxation test, the sample was assumed to 
ie be in simple shear. The shear strain (7%) was calculated as the tangent of 
ie Reference 2 - Jones, J., "Solid Propellant Structural Integrity Investigations: 
a4 Dynamic Response and Failure Mechanisms in Solid Propellants," RPL-TDR-64-32, 

Vol. I, Lockheed Propulsion Co., February 1964. 
: 
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the displacement angle or AL/G.L. The shear stress (T) was calculated as 
force/area (area = 3 sq. in.). 


The data modification to insert peak and minimum stress points previously 
discussed were utilized for the biaxial and nonisothermal tests. 


3.2.1 Biaxial Constant Rate Test No. 14 

The biaxial constant rate tests to failure were run with the 1/4 x 1-1/4 x 
6-in. specimens of UTP-3001 and UTP-19,360B. The 40, 70, and 120°F tests were 
run at crosshead rates of 2, 0.2, and 0.2 in./min. The test equipment, with 
specimens in place, is shown in Figure 49. A typical load-time curve is shown 
in Figure 50 for UTP-19,360B at 71°F and 2 in./min. crosshead rate. Because of 
the fixtures and more difficulty in adjusting linkage than with the 6-in. bar 
specimens, the three samples did not start loading simultaneously. Sample 2 was 
adjusted to an effective zero and is shown in Figure 51. Tabular data are given 
in Table 24. 
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Figure 41. 1/2-in. Bar Stress Relaxation Data at 3% Nominal Strain 
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3.2.2 Biaxial Straining-Cooling Test No. 15 

Biaxial specimens of UTP-3001 and UTP-19,360B propellants were simultane- 
ously strained and cooled from 115 to 40°F at a crosshead rate of approximately 
3.x 107> in./min. over a 40-hr period. The results for UTP-3001 are shown in 
Figure 52. The stress-time traces for all three samples appear to start 


together but spread out as the test progresses. The tabular data are given 
in Table 25. 


3.2.3 Biaxial Stress Relaxation Test No. 16 

Biaxial stress relaxation tests were run with the 1/4 x 1-1/4 x 6-in. 
specimens of UTP-3001 and UTP-19,360B propellants at a nominal 3% strain and tem- 
peratures of 40, 70, and 120 F. The 4o°F data for UTP-19,360B are shown in 
Figure 53 as typical with good reproducibility. The loading ramp rate was 
0.2 in./min. Data are given in Table 26. 
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Biaxial Biaxial samples of UTP-3001 and UTP-19.360B were ramp loaded to 
constant rate failure at rates of 2. 0.2. and 0.02 in./min. at temperatures of 
41. 70. and 120°F. 


Biaxial Biaxial samples of UTP-3001 and UTP-19,.360B were simultaneously 
straining- Strain and cooled from 120 to 40°F over a 40 hr. period. 
cooling 


Biaxial Biaxial samples of UTP-3001 and UTP-19.360B were run in stress 
relaxation relaxation tests at 40. 70, and 120°F. 
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Shear Shear samples of UTP-3001 and UTP-19.360B were run in stress 
relaxation relaxation tests at 70°F. 


6-in. bar 6-in. bars of UTP-3001 and UTP-19.360B were simultaneously strain 
Straining- cooling and cooled trom 120 to 40°F at three slow rates 
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Biaxial Biaxial samples of UTP-3001 and UTP-19.360B were run in ramp- 
thermal relaxation-ramp tests with simultaneous cooling or heating {i.e.. for 
similitude reverse ramp) 


Note: Nominally three samples were run for each test and condition. Legend: 
T = temperature 
t = time 
€ = strain 


Figure 42. Biaxial and Nonisothermal Phase III Testing 
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3.2.4 Shear Relaxation Test No. 17 
Shear relaxation tests were run on 1 x 1 x 3-in. samples of UTP-3001 and 
UTP-19,360B propellants. The samples were post bonded to steel plates and run 
one at a time by loading them at 0.2 in./min and ambient temperature with offset 
fixtures so the load was transmitted through the center line of the sample. The 
3 samples for each propellant were hand reduced and digitized for computer stor- 
age and printout. A typical example is shown for UTP-3001 in Figure 54 with 
a data given in Table 27. Peak stresses and strain were very close as was the 
1 hour relaxation stress on each propellant even though the samples were run 
separately. 
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is! 3.2.5 Straining-Cooling Multiple Rates Test No. 18 
The rate effect on the straining-cooling response was determined on 
i UTP-3001 and UTP-19,360B propellants. The 1/2 x 1/2 x 6-in. bar sample was used 
os so that testing could be completed in the shortest time possible. The rate 
ay effect for the uniaxial specimens was then applied to the biaxial test No. 15. 
Cooling was from 110°F to 40°F at the crosshead rates of 0.002, 0.0002, and 
“ 0.0004 in./min. Typical data are shown for the 0.002 in./min. rate for 
UTP=-3001 in Figure 55. Good reproducibility is shown within the set of 
3 samples. Data are given in Table 28. 
a 3.2.6 Biaxial Quinlan Complex History Test No. 19 
on The 1/4 x 1-1/4 x 6-in. biaxial samples of UTP-3001 and UTP-19,360B were 
ie subjected to the complex cycling and relaxation history indicated in Figure 42. 
When the tests were run the linkage, misalignment, etc. was such that the 
5 samples were not loaded an equivalent amount. The first sample to be loaded had 
s the correct strain determined but strain on the other samples had to be adjusted 
fa to the time the stress ramp started on each. Each sample was reduced separately 
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Figure 44. Principal Stress and Strain Ratios at the Center of Biaxial Strips 
for Varying Height-to-Width Ratios for a Poisson's Ratio of 1/2 aa 
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and data were modified to pick up the peak and minimum stress points. Sample 1 
for UTP-3001 is shown in Figures 56 through 58 where the complex test has been 
divided into segments on an expanded time scale to show test details. The first 
cycle in Figure 56 showed no load and the second cycle showed very little. By 
contrast sample 3 (not included here) had a first peak stress-strain of 37 psi, 
2.26% and second peak of 76 psi, 4.97%. During the unload part of the cycle 
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Figure 45. Normalized Axial and Lateral Strains Along the Midplane 
Biaxial Strip Specimen 
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after stress reached zero, strain decay was estimated from other tests where 
strain was measured by cathetometer. Data are given in Table 29. 


3.2.7 Cyclic Testing Test No. 20 

The complex cyclic testing of UTP-3001 and UTP-19360B propellants was done 
with the 1/2 x 1/2 x 6-in. bar samples. These tests were run at 70% with 
cycling from 8 to 4% then 12 to 8% followed by 8 to 4% strain at a crosshead 
rate of 0.1 in./min. The strain levels were set so that the sample would 
not reach zero stress on the unload cycles. By doing that a correct evaluation 
of strain was obtained during the tests. Data for the UTP-19,360B ri=apellant 
are shown as typical in Figure 59 with digitized results given in )...1# 30. 
This propellant broke without completing the nominal 30 cycles but UTP-3001 
went all the way through the test. The nominal 10 cycles per segment was 
dependent upon when the time scale matched available personnel. It was 11 
cycles for the first 2 segments for UTP-19360B as shown in Figure 59. 
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Figure 46. Strain Variations Along Midline of Strip Specimen 
(N = 0) for Poisson's Ratio of 1/2 
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3.2.8 Biaxial Ramp-Relax-Ramp Test No. 21 

The ramp-relax-ramp tests were run on UTP-3001 and UTP-19,360B propellants 
with the 1/4 x 1-1/4 x 6-in. biaxial specimens. The first test was ramped at 
0.0005 in./min. to 6% strain and simultaneously cooled from 120 to 70°F. It was 
held at 6% strain nearly 23 hours then ramp to failure while cooling towards 
4o°F. The data for UTP-3001 are shown in Figure 60 with data points given in 
Table 31 as shown in Figure 60. The cooling cycle did not end at the peak 
strain consequently the relaxation of stress was not the normal type behavior. 
The continued cooling increased the propellant stress capability so that 


the normal relaxation behavior did not start until the propellant temperature 
stabilized. 


This test was repeated starting at 110°F and taken to 6% strain with a peak 
stress of 70 psi compared to 30 psi for the above test. The longer ramp time 


120 


0.6 


a, ando,/o, 


a 


a 
° 
> 


ae 0.2 


& = x/a 
Figure 47. Stress Variations Along Midline of Strip Specimen 
(N = 0) for Poisson's Ratio of 1/2 
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im 

“ allowed the cooling to reach 40°F at the peak stress. The samples were allowed 

5 to relax overnight and then unloaded to 3% while warming the samples to room tem- 
a perature. Figure 60 was considered to be sufficient to represent the test. 

: 3.2.9 Propellant Aging Effects During Phase III Testing 

am The biaxial testing in Phase III was scheduled to minimize the dry box stor- 
" age time after sample machining. The purpose was to reduce the within carton or 
3 box variability encountered during the wiaxial testing in Phase II. The major- 

re ity of the testing was accomplished within 1 week of machining the samples. 

fae Those that went over a week did not seem to be influenced excessively. 

a 

. The data obtained for initial ramp loadings (i.e., undamaged behavior) on 

— the tests run are given in Table 32 for UTP-3001 and Table 33 for UTP-19,360B. 


oe The biaxial constant rate data for UTP-3001 are plotted in Figure 61. These 
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Figure 49. Biaxial Test Setup and Instrumentation 
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results were used to obtain comparison data for other tests at different strain 
levels in Table 32. The stress at 8% strain versus temperature is shown in 
Figure 62 for UTP-3001. The tests at different temperatures were taken from 
different redwood boxes of propellant and there appears to be some between box 
and sample differences. The data shifted for rate effects to the 0.2 in./min. 
in the lower plot indicates that extrapolation of the temperature stress plot 
would be unreliable. This eliminated any reasonable direct comparison with the 
Slower rate straining - cooling tests as noted in Table 32. The comparisons did 
show reasonable agreement with some samples from each of the redwood boxes used. 
The straining-cooling data are expected to be of the same general family. 


(Text continued on page 168.) 
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Test No. 17 - Stress While Step Straining for UTP-3001-750/7768 


Figure 54. 
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Figure 61. Biaxial Constant Rate Data for UTP-3001 
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The UTP~19,360B biaxial constant rate data plotted in Figure 63 were used 

to obtain the stress comparison data (Table 33) at different strain levels. 

The stress-temperature plot in Figure 64 gave a good temperature extrapolation 
which was used to obtain comparison stress values for the slower rate straining- 
cooling tests. The stress values selected for comparison were early in the 


test at an elevated temperature and small strain. This was chosen to specifically 


avoid the complication of evaluating the combined thermal-mechanical interaction 
shift factor (Ap). While one of the straining-cooling tests showed a difference 
of 35%, the absolute delta stress was small. 


Uniaxial, Biaxial and Shear Comparison 

Comparisons of the different samples were made in order to show that the 
propellant, used in each of the tests, was of the same family. This was done at 
ambient temperature for selected rates and was limited to strain rates that were 
close to each other. This minimized the time-temperature equivalence shifts to 
small changes for neglectable data input errors. The adjustments made for 
strain levels are given in Table 34. 


The comparisons between uniaxial and biaxial in the table are close to the 
theoretical ratio of 75 to 80%. The shear to uniaxial ratios of 0.28 and 0.39 
bracket the nominal theoretical value of 1/3. For the comparisons made in Table 
34, the UTP~3001 and UTP-19,360B propellants have to be considered part of the 
same family. Any minor differences can be attributed to a between carton 
effect. 
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Figure 63. Biaxial Constant Rate Data for UTP-19,360B 
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4.0 THEORETICAL DEVELOPMENT 


4.1 INTRODUCTION AND PRELIMINARY STUDY 

Generally speaking, solid propellants may be considered as lightly cross- 
linked long-chain polymers, highly filled with coarse solid particles. They 
respond viscoelastically to the action of external stimuli, but certain aspects 
of their behavior cannot be reproduced by the classical linear or nonlinear 
theories of fading-memory materials. Thus, in recent years, much work has been 
concerned with the development of appropriate models to predict the mechanical 
response of solid propellants. A current trend is to express the observed 
response in terms of some measure of "damage" at the continuum level where 
damage is described as the difference between the observed response and that 
predicted by a fading-memory theory, like Linear Viscoelasticity. There is now 
sufficient experimental evidence to show that damage (References 13, 17, 20, 24, 
and 28) per se is a microscopic phenomenon associated with the initiation and 
growth of flaws, debonding between matrix and solid filler particles, and molecu- 
lar chain scission. Although it is largely irreversible, damage is partially 
recoverable shortly after removal of the loading system. This recovery from 
damage is termed healing. It is clear that, depending on the propellant and 
service requirements, it may also have to be accounted for in a constitutive 
theory for solid propellants. 


In the present program, two approaches to characterizing damage have been 
followed. In the first one, damage is treated as the algebraic difference 
between the measured stress and that predicted by Linear Viscoelasticity, so 
that: 


Go (t) = O(t) - og (t) (1) 


in which t) and Og are the linear-viscoelastie and correction terms, respec- 
tively, witho, the measured stress. In the second approach, the difference 
between measured and fading-memory type stresses is handled through a stress- 
correction function in the following form. 
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(2) 
a(t) s Cc (€ maxs eee) o g(t) 


The softening function (C) is made to depend on the past maximum strain or 
stress and 0 p(t) represents an appropriate function of the fading-memory type 
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Broadly speaking, the current versions of models by M. Gurtin and 
M. Quinlan are of the type presented in equation (1) above, while those of 
R. Schapery, W. Hufferd and Swanson are of the form given by equation (2). 


The following presents some experimental evidence on the type of non- 
linearities exhibited by solid propellants, and briefly discusses the pioneering 
work of Mullins and Tobin in treating the large hysteresis observed in tire 
rubbers. Next, the theory of Linear Viscoelasticity is applied to predict the 
response of UTP-19,360B and UTP-3001 under various strain histories. The 
ensuing results are meant as a basis for comparing the propellant response as 
predicted by each of the candidate constitutive laws. This comparison should be 
most meaningful because each of the theories considered evolved from a set of 
modifications to Linear Viscoelasticity. Subsequently, the nonlinear theory of 
Farris (Reference 5) is presented. This theory was employed during the first 
phase of the program to predict the response of TP-H1011 and to compare the 
results with those of the other five constitutive laws. Finally, a detailed 
description is given of each of these five stress-strain relations. This 
includes the original concept of the models, their current versions, comparisons 
of predicted and measured stresses for a variety of strain histories, and some 
pertinent guidelines for characterizing solid propellant according to each 
theory. 


4.1.1 Experimental Background 

The complex behavior of solid propellants, as well as some attempts at 
developing usable stress-strain laws for these m-‘erials, are well documented in 
References 3 through 31. It is shown there that a given deformation process 
causes a change in the response properties of solid propellants, for instance a 
drop in the relaxation modulus. As stated before, this deviation from some 
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expected response is what has been called damage. It is evidenced as 
phenomenological macroscopic changes that are caused by undefined, but real, 
irreversible or partially reversible microscopic changes. Polymer bond 
breakage, vacuole formation in the polymer matrix, dewetting between the polymer 
matrix and solid filler particles, microcracking, etc., are among the possible 
microscopic causes of observed permanent-memory effects in propellants. 


Studies on uniaxial solid propellant samples have indicated that these mate- 
rials exhibit large hysteresis even at small strains. These studies have also 
revealed that the state of damage in solid propellants is determined primarily 


by the maximum strain or stress undergone during the loading histories. 


The typical nonlinear hysteresis and permanent-memory effects exhibited by 
solid propellants are illustrated in Figure 65. A series of finite-duration, 
variable-strain-level ramp pulses were used to obtain the propellant response 
subsequent to a given damage history (References 12 and 13). All ramps had the 
same initial moderate rise rate, with two exceptions to be noted later, and all 
ramps had the same very slow decline rate. 


Observations of the load on the specimen after returning to its original 


length (zero strain) showed that it took approximately 30 min for the stress to 
relax to zero. 


A series of tests were run on a 1/4 x 1/4 x 4-in. tab-end sample. The 
virgin specimen was initially strained to a level of 7.04% and allowed to relax 
to achieve a rest-state condition. The first part of the testing is shown in 
the lower half of Figure 65 (curves A-H) and the last part in the upper half 
(curves H-M). 


Curve A shows the load response to this first pulse. The specimen was then 
subjected to four successive ramp strain pulses ranging from €° = 3.82% to e° = 
6.34%. There was a rest period allowed between each pulse. The results are 
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shown in Figure 65 as curves B through E. 
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Figure 65. Relaxation after Damage 
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Two aspects of the propellant's behavior are worth noting. First, during 
the constant strain rate portion of the ramp, each successive load-time curve is 
essentially identical. This indicates that the "new material" has the same 
nonlinear rate-dependency under repeated strain conditions as long as the strain 
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levels are below the previous maximum strain of €° = 7.04%. Second, the relaxa- 
tion portions of the curves are essentially homologous, indicating that a 
viscoelastic relaxation process is taking place. 


Curves F and G present the results of two additional tests at two succes- 
sively lower strain rates where the sample was loaded to 5.64% strain each time. 
A strong rate dependency is observed during the rise portion of the ramp; 
however, curves F and G rapidly rejoin curve D indicating that the material is 
behaving in a viscoelastic fading-memory fashion. 


The specimen was next subjected to a ramp strain pulse reaching a higher 
strain level (€° = 9.86%) than the maximum 7.04% previously experienced (Fig- 
ure 65, curve H). The first part of curve H repeats the loading ramp portion of 
curves B-E to indicate the same "new material" rest-state. Note that the load- 
time curve returns to the initial or virgin constant strain rate curve once the 
previous maximum strain (7.04%) has been passed. 


Subsequently, the specimen was strained with the ramp pulse to four 
different strain levels less than 9.86% (€° = 4.26%, 5.64%, 7.04%, and 9.01%), 
as shown in curves I through L. The results show that a new rate~dependency has 
developed (compare the rising portions of I through L with the rising portion of 
H). Thus, another "new-material" rest-state has been produced as a result of 
the second maximum strain level of 9.86%. Lastly, the specimen was strained to 
another new maximum of €° = 12.68% as shown in curve M. It again appears that 
it returned to the virgin undamaged curve once the 9.86% strain level was 
exceeded. 


The above experimental evidence suggests that the form of the constitutive 
equation should remain unchanged with respect to the material's current rest- 
state. This condition should remain as long as the damage is unchanged (i.e., 
the € max is unchanged during its subsequent strain histories). 


Figure 66 shows a replot (curves N and 0) of some of the results just 
discussed. After an initial maximum strain (7.04%) the specimen was allowed to 
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Figure 66. Relaxation after Damage 
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re return to a rest~state and then strained to a value of €° = 5.64%, with the 
ane result shown as curve P. These three identical strain history tests of three 
ae different material states indicate that the higher the state of damage 
(primarily €max), the softer the material response upon subsequent testing. 


er In addition, other experimental studies have pointed out the importance of 
of healing effects, load duration, and initial strain rate (References 12, 13, 14, 
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and 15). Finally, it is important to note that the behavior of solid propel- 
lant, depicted in Figures 65 and 66, is similar to that exhibited by rubber. The 
nonlinear uniaxial stress response of rubber, with and without carbon black fil- 
ler and in the absence of time effects, was characterized quite well by Mullins 
and Tobin (Reference 27) with equation (3). 


where: 


€= engineering strain. The Mullins-Tobin model is not limited to small 

strains. 

€y =€y(o) = strain as a function of engineering stress for the polymer without 
filler and without damage. The characteristic shape of this function is 
shown in Figure 67. 

F = F(o,,,, N) = damage or softening function which depends on the maximum 
stress experienced by the rubber and the number N of loading and 
unloading cycles. F is not very sensitive to N, but depends strongly on 
any hard filler particles that may be present. 


A large amount of rubber data can be predicted by means of this equation 
when the samples are not allowed to rest between cycles. Recovery or healing 
occurs as a function of the rest time. Therefore, healing would have to be 
considered in an accurate characterization of rubber. 


Introducing the inverse of €,, = €,,(0), equation (3) may be put in the form: 


€ =f (€/F) (4) 


which shows that F (where F€1) is a strain-magnification factor. The ratio €/F 
is interpreted by Mullins and Tobin to ‘be the average strain in the rubber phase 
of a hard particle-filled rubber. Without damage in a highly-filled rubber, 
F<<1. As the rubber is cycled between the strains €= 0 andeés€may, the ratio €/F 
at any strain decreases, and therefore the stress decreases. The shape of the 


stress-strain curve is still as shown 
in Figure 67. It is similar to that 
for solid propellant after first-time 
loading. This fact and the ability of 
the model presented in equation (3) to 
reproduce a large amount of rubber data 
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explains the great influence of the 

Mullins-Tobin approach on the devel- 
Strain opment of nonlinear constitutive 
theories for solid propellants. 
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Figure 67. Stress-Strain 
Curve for Rubber 
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4.2.1 Linear Viscoelastic Constitutive Equation 
4.2.1.1 Linear Viscoelastic Model 
The one-dimensional stress-strain law for a thermorheologically simple 
linear viscoelastic solid may be expressed as: 


t 
dt 
o(t) = | E (St - S-) ——(1) dr 


dt 
fe) 


stress 
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strain 
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relaxation modulus (PRONY series representation 
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using a matrix solution for curve fitting data; 
CSD Data Analysis Procedure No. 7.3) 
= temperature-shifted time, given by: 


as 
ay [T (7) ] 


ar = time-temperature shift function, taken in the form: 


ay ae 


on 


yyy 4% 


Paige 


a 

» Tr - Ta (7) 

= Ap = | ———— 

bi T-Ta 

es in which Tp is the shift reference temperature, and both Tg and m are material 

ne parameters. The material parameters were obtained using CSD Data Analysis 

t Procedure No. 7.4, which is a curve fit routine using Powell's algorithm. 

x The linear viscoelastic model was used to predict the response of 

. UTP-19,360B and UTP-3001 under several strain histories. The corresponding 

= results are included here as a basis against which to compare the stress predic- 

. tions obtained using the nonlinear stress-strain laws considered in the program. 

x! 4.2.1.2 Stress Predictions 

a The measured response is compared against that predicted by linear 

viscoelasticity for UTP-19,260B in the following order (Figures 68 through 78). 

os The results for the lowest and highest constant-rate tests (Test No. 1) 

a appear in Figures 68 and 69. Those for the dual-rate tests (high-to-low and 
low-to-high, Test No. 3) are shown in Figures 70 and 71. Figure 72 contains the 

a comparisons for a saw-tooth strain history (Test No. 5) with increasing strain 

cs peaks and rest periods between cycles. The results corresponding to short- and 

SS long-duration similitude tests are presented in Figures 73 and 74. The pre- 
dicted and measured responses for the three-step relaxation test (Test No. 13) 

x are included in Figure 75. In addition, the time-temperature superposition prin- 

s ciple is put to use with constant rate tests (Test No. 1) at 70, 40 and 123°F, 

yy as shown in Figures 76 to 78, respectively. Stress predictions for some of the 

a same tests on UTP-3001 are shown in Figures 79 through 86. 

ES: 

» 4.2.2 R. Farris Nonlinear Theory for Solid Propellants 

: i The work of R. Farris (Reference 5) was a major attempt at predicting the 

_ nonlinear response of solid propellants in rocket motor analyses. 

¥ , Experience with Farris’ stress-strain law at Chemical Systems showed that 

Re this theory could not predict the response of solid propellants under strain 

* (Text continued on page 200.) 
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histories that were not included in the set used for material characterization. 
Although the predictions were generally acceptable for loading histories of the 
types included in such sets. 


The nonlinear theory of Farris was considered in the first phase of the 
program and it was compared to the other five approaches originally proposed. 


Farris postulated a model, based upon previous work on rubber elasticity, 
to account for the permanent-memory effects exhibited by many solid propellants 
under constant rate loading. This model presumes the existence of inhomo- 
geneities in the local strain field between filler particles, a distribution of 
polymer chain lengths between filler particles, and a uniform failure strain for 
each polymer chain. This model has been successful in predicting the nonlinear 
permanent memory response of solid propellants before dewetting. The model 
predicts the same response in compression as in tension. This prediction is in 
agreement with experimental observations, although the molecular mechanisms con- 


tributing to the permanent memory response in compression are clearly different 
from those in tension. 


Once dewetting occurs, the model is modified to account for vacuole forma- 
tion and different results in compression and tension are expected. Farris 
presented the constitutive equation as the sum of an essentially time- 


d 
independent bulk stress op and a time-dependent deviatoric stress 0; j» so that 
in general 


d 
o1j (t) sop 545 +053 (t) (8) 


where 5, j is the Kronecker delta equal to wmity if i = j and zero otherwise. The 
form developed for the deviatoric stress is 


m2 
d BI d Ty d 
043 (t) = e-¥id/Iy Ay eqs (t) + Ag 7? eqy (t) (9) 
‘Y"P2 
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.d Ty my d 
: f Ay (t = &) ayy (E) ab + Ay (t= €) B45 (LD a 
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where 
Ig = volume dilatation = e141 + e292 + €33 for small strains 
1 2 2 10 
I, = octahedral shear strain =—[len - e202) + (e092 = €33) ote 
3 2541/2 
+ (e33 - e414) ] 
d 
eyj = deviatoric strain tensor 
B, Ay, Ao, Ag, Ay, Mo, My, Po, Py, = constants 
and 
t 1/P 4 
Py 
IIyIlp, = flac! dt (11) 


° 


Equation (9) has been applied to reasonably complex deformation histories using 
unpressurized and pressurized uniaxial and biaxial test specimens. The 
agreement was not as good as would have been desirable, but it was still better 
than Linear Viscoelasticity. Time-temperature superposition was included in 
equation (9) by introducing a time-temperature shift factor, at, and redefining 
the Lp norm of equation 140 in the form: 


t 7 


WIyll,, = ag 
‘ Pi J ar(&) 
° 


Experimental data for simultaneous cooling and straining have been fit using 
equation (9) with the introduction of a time-temperature shift function ay 
through equation (11). The justification for introducing an ar in the above 
manner is not immediately obvious or adequately explained in the available lit- 
erature. 


To present the response to interrupted and cyclic constant strain rate 
tests, equation (9) was modified by setting Py = oo and Az = -Ay 90 that: 
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d { ly m2 4 
EN (t) = e ~Blyg/Ty | Ay eqs (t) + Ag (——— eyj (t) (13) 
; NIy lp, 


Iy \"4 a 
+] 1- Az (t - &) eqs (&) dé 
ITy Woo 


I Lyllog# nax|t 4(8)| » O< Ect 


The multiplier for the hereditary integral in equation (12) vanishes 
whenever the current value of I, is at its largest, and is non-zero for all 
other values. This representation allows for viscoelastic (fading memory) 

response on unloading. 


The bulk stress, 0, in equation (8) was taken to be essentially time- 
independent, although there is coupling between distortion and dilatation as 
indicated in the exponential multiplier in equations (9) and (13). 


The first attempt to represent the bulk stress took the form of a series 


o 


(14) 


N ‘oj 
’ Crk = EY = Aggy Tg Ty § Ao 
j20 


This equation adequately predicts the bulk stress as long as it does not vary 
greatly. However, when a hydrostatic pressure is superimposed, very poor 
r results are obtained. In an attempt to overcome this difficulty, Farris modeled 
the compressibility of the gas voids caused by vacuole dilatation by treating 
them as spherical voids contained in an elastic medium. Assuming that (1) the 
voids themselves offer no resistance, (2) that all void dilatation is caused by 
distortion of the surrounding elastic material, (3) the void content at zero 
pressure may be represented as a power law in terms of the octahedral shear 
strain Iy , and (4) that the bulk behavior varies linearly with hydrostatic 
pressure (P), the model yields: 


ee ae oe 


wet ca BASES Ved Sogn hen Sy PR Rag eat at a RLS Som i Ck Siig? Sag oO Sd en ge 
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for the dilatation. The C;, Co and n are constants and G is the shear modulus 
of the elastic matrix material. 


ee, 
TIPTREE 


4.2.3 R. Schapery's Nonlinear Stress-Strain Law 
4.2.3.1 Original Model 

The constitutive theory advanced in References 17 and 18 for viscoelastic 
materials with microcracking was taken by Dr. R. Schapery as the starting point 
for predicting the response of solid propellants under general loading condi- 
tions. The one-dimensional version of this law takes the following simple form: 


ane 


Bete tee 


Sirageg 


(16) 


where op is the linear viscoelastic stress for a thermorheologically simple mate- 
rial: 


d 
fee-mter, 
dr 


° 
= € - a(T - Ty) = strain due to mechanically applied 


stress 
t 


| at'/ag [T(t") ] 


° 


ft = & (7) 
E(&) = linear viscoelastic relaxation modulus 


To = temperature at t = 0 


ae 


age 
ar 


Ap = Ap (T) = temperature-dependent material function 


~@ 
Jove fe: 


A= \ (Sg): softening function in which the damage parameter: 


See 9K 


Se = 


depends only on the strain and temperature histories, and: 


1 for O€ € <€] 
(e/ey)* for €1<€ <9 


(eo/er)? for € > €2 


for constant threshold strains €; and €5, withf>0. 


The function F = F(€g) and the positive, constant exponent, q, originate 
with the equation for microcrack speed, 


dA q 
— is M(Ky/f) 
af 


where M is a positive constant and: 
df = dt/Ag 
in which A, = A,(T) is the shift-factor for microcrack growth rate. 


The functional form of the softening function, 4 = A(Sg), depencs on the 
type of behavior that need be reproduced. The following special case was used: 


refi + eSo p/q 


where c and p are positive constants. Note that when Sp = 0, or c = 0, a linear 
viscoelastic stress-strain equation is recovered from Equation (16). 


Taking Ap = 1, several sets of numerical values for the constitutive parame- 
ters corresponding to TP-H1011 were tried without success. This theory was also 


used to predict the response of UTP-19,360 and UTP-3001. Having failed to 
perform better than Linear Viscoelasticity in many cases, it has undergone sev- 
eral changes since. 


4.2.3.2 Current Model 
The essential form of the modified uniaxial stress-strain relation is given 


by: 


o = f (€°, €°m, S) 


where: 
= engineering stress 


°: psuedo strain 


t 


1 de 
=o—_ E (t - 7)——dr 
ER dar 


° 
= maximum value of |€°/ up to the current time 


= damage parameter 
t 
q 
= f |e>| dt 
° 


ER = arbitrarily selected reference modulus, 


1/q 


E(t) = linear viscoelastic relaxation modulus, 
= Eg + Ep t"? & Ep (E, + t7"), 
Er = Ee/Eo, 


q = positive constant. 


The functional form of f in equation (21) depends on the material considered. 
Studies on solid propellant to date indicate it may be taken as follows for some 
solid propellants. 


f = Y; Yo Y3 Py5 sign (€°) (25) 


in which 
1 if 6° >0 
sign (€°) = Oo if €9=0 


-1 if «9<0 


and Pis is used to normalize function Y3 to unity, at a reference point. The 
Y,"s are the following functions of damage and pseudo strain: 


1 + AyS + ApS? + A383 for S < Sp 
Y, = Y, (S) = (26) 
Ay S45 for S >S, 


(27) 
Yo = Ag $0+63-Sx (em) (0-463 - My - Ly) |e? ie 


¥3 = Co + C4x + Cox? + Cgx3 + Cyx +C5x5 (28) 


where: 


(29) 


in which x, 1s the only root of the equation: 
max(S,) = Y3 (xp) (30) 


with max (S,) representing the maximum value of S,, up to the current time, and: 


206 


-~e+n? 


SREB RE PR San Aa Tees Sa PL SE Be AB COP PM SG Te oe Dea MATa: tee Ee ante SES Nie, Sai eMC oes aie ae Stay wa en Be 


(31) 


while A is a factor that accounts for relatively small higher order effects 
possibly due to rehealing and particle interaction. 


The constants entering the definitions of Y+, Yo, and ¥3 depend on the 
material. For UTP-19,360B they are: 


42 (32) 
0.637 

-0.387 

0.85 


and the factor Ais given by: 


a 


ieke ca (33) 


The resulting form of equation (21) for UTP-19,360B is thus: 


L 


x 
o= Py5 Ag Yy Y3 |< Sign (€°) (34) 


Clearly, if Ly = 1, equation (34) may be written as: 
t 


8 he [fes0— 3 (35) 
aT 


~~ wra ee 


° 
in which ap = ap (e°, €9, 8) plays the role of a softening function, reminiscent 


of the Mullins-Tobin approach. 


4.2.3.3 Stress Predictions 


The current version of the nonlinear model developed by R. Schapery may be 


used to predict the response of solid propellants with a rather remarkable 


degree of accuracy, as may be seen in Figures 87 to 94, which are sample cases 


of the isothermal tests considered in the program. The first two plots 
(Figures 87 and 88) correspond to the highest and lowest constant rate tests 


(Text continued on page 216.) 
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(Test No. 1) available in the data for which the difference between theory and 
test is greatest. Figure 89 shows the saw tooth test (Test No. 5) at constant 
rate with increasing strain peaks. Figures 90 and 91 pertain to the dual-rate 
tests (Test No. 3) Results for the short- and long-duration similitude tests 
(Test No. 12) are given in Figures 92 and 93, and Figure 94 includes a 
three-step relaxation test. 
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“ Finally, it is important to mention that a complete characterization of 
UTP-19,360B was also carried out using Ly = 1 (the value leading to equation 
(35)), and the ensuing response predictions were very close to those obtained 
with L, = 0.85; only the low-to-high dual rate test of Figure 100 was predicted 


4 

ey somewhat better with L, = 0.85. 

) 

“4 

4.2.3.4 Material Characterization 

ee In evaluating the material constants and property functions, the following 
a observations may be valuable: 

id 

os (1) Y, appears to be a decreasing and concave down function of damage, as 
“| presented in the following figure, its variation being brought about by 
ed 


vacuole formation. 


ores 
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(2) The function S,, which provides a certain measure of damage, increases 
as a direct result of a reduction in the number of polymer chains 
supporting the internal stresses; the larger S,, the higher the stress 

ms! on each chain. 

rs (3) For UTP-19,360B, the state of damage is essentially constant during 

E unloading and reloading, and the shape of the so-called backbone curve 

resembles the stress-strain curve for rubber, which is of the form: 


ere 
eS otae sees 


DO ee TE SL ST i i a ae Pr Sg oe RE gS ee Tg eee en Ea SCRE S Le OL Vat ees | eee ie! 


=, 

= 

a 0.85 

a x3 |€0| (36) 
f 

xe as shown in the sketch below, in which the steepness increased with 
‘Te increasing Sy. 

~ 0 

& (1) Loading portion 

aa (2) Unloading and reloading portion 
ay (backbone curve) 

: “6 

ta 

oo (4) The selection of L, can be made by plotting unloading data in the form 
na suggested in the following diagram. 

= —~— Unloading from point 

om on virgin curve 

a Ly =0.85 

a €) Lx 

2 ee 

of Noting that the quantity: 

“ g 

by (eo) bx 


resembles a secant modulus, and that for most tests of UTP-19,360B 
eo. Ly = 0.85 produced a finite limiting value as €° approached zero, it is 
x Suggested that L, be found in this fashion for other propellants. 


se (5) For constant-rate tests, one has: 
os 

C eo =? 
~. m 
ae N = 1 
i 

ae and thus, from equation (29): 

pr X= X, 
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Also: 
Sp = max (S,) 
from which: 
Y¥3 = Sp 


Equation (34) then reduces to: 


o = 1.861 Yy g0 .637 (€ 00-463 


(37) 
with Y; = Y; (S) given by equation (26). 
For very small damage: 
Y, (S) = 1 
so that equation (36) becomes: 
o = 1,861 $0-637 (¢0)0-463 (38) 


in which the stress increases with damage, probably because of molecular 
chain stiffening due to an increase in stress per chain. 


With the foregoing observations in mind, determining the material 
properties can be accomplished as follows: 


(1) The exponent, n, appearing in the relaxation function, is obtained 
from relaxation-modulus data. 

(2) The normalized coe:ficient, E,, entering the relaxation modulus, 
is determined to make unloading curve 2 in the figure above pass 
through the origin. 
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(3) The exponent, q, present in the definition of the damage parameter, is 
evaluated using equation (37) and two constant-rate tests at small 
values of damage. 

(4) The function Y; is obtained by curve-fitting equation (36) to 
constant-rate tests over all strains out to failure. 

(5) Experience to date indicates that the function Yo is independent 
of S and €,, and therefore equation (30) may be used instead of 
the more general form of equation (37). 

(6) The backbone curve Y3 is determined using wmloading and reloading data 
like those available in a cyclic test whose first peak strain is the 
largest. 

(7) Finally the correction factor, ’, can be ascertained from a relaxation 
test at a large strain level. 


4.2.3.5 Multiaxial Generalization 

A micromechanics model has been developed which predicts the form of 
equation (30), and it is presently being used to develop a multiaxial form 
of the theory. 


4.2.4 M. Gurtin's Theories for Nonlinear Viscoelastic Materials 

Four essentially different approaches have been followed by M. Gurtin in 
trying to predict the response of solid propellants that exhibit damage. The 
stress-softening theory appears to be the most accurate of the four laws as will 
be pointed out. 


4.2.4.1 Original Model 

The one-dimensional stress-strain law for materials undergoing internal 
damage was based on the hypothesis that the state of damage at any time is 
completely characterized by the maximum strain, e,, that the material has 
experienced. 


€m (t) = max €(s) (39) 


Onas et 


ac pe por, 
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The stress, 0, is given by a constitutive equation of the form: 


a(t) = gle(t),en(t)] (40) 


and it depends only on the current values of strain and damage. Such an 
equation is, of course, rate-independent. 


In this theory, if the maximum strain occurs at the present time, then: 


Em (t) = e(t), (41) 


and equation (39) reduces to: 


o = Gem) = elem €m) (42) 


The stress-strain curve: 


is called the virgin curve and is traced out in an experiment with monotonically 


increasing strain. 


Using the virgin curve, equation (39) may be rewritten in the form: 


o = F(E,€m) Gem) (44) 


with: 


ae Be (45) 


the relative strain, and: 


g( Em, € ) 
ti. = (46) 
G(E m) 
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The function F(§, em) is called the damage curve at the damage level €m, and is 
such that: 


F(1, Ep) = 1 (47) 
In some situations of interest F(£, €,) is independent ofe »: 
F(E,€m) = FC(E) (48) 


when this is so, F(&) is referred to as the master damage curve, and equa- 
tion (43) reduces to: 


o= F(E)G(Em) (49) 


As pointed out previously, this is a rate-independent theory, and as 
such, cannot be used for loading rates that differ much from that used to 
determine the damage function. This situation was remedied by changing the 
stress-strain law to the one described next. 


4.2.4.2 Nonlinear Model Based on Stress Softening 
To develop a simple theory of stress softening which allows for rate 

effects and which returns to Mullins’ original idea of using the past stress 
maximum as the damage parameter, two fundamental ingredients are considered. The 
first is the virgin stress, S, which represents the stress the material would 
experience in the absence of softening. This stress is assumed governed by a 
constitutive equation of the type encountered in linear viscoelasticity. The 
second ingredient is a damage function, F, which gives the true stress, 0, when 


the virgin stress, S, and its past maximum S are known. 


The one-dimensional form of the constitutive law for a classical linear 
viscoelastic material is given by: t 


o (t) = fow-nien dr 


—oo 
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in which o (t) is the stress; ¢(t), the strain; and G(t), the relaxation 
function. It is further assumed that ¢(t)=0, prior to t = 0. 


The generalization of equation (50) is begun by defining the quantity: 
t 
S(t) = foe - 7) €(r) dr 


=-T 


(51) 


which is called the virgin stress and which represents the stress that would be 
present in the absence of softening. It is assumed that the extent of softening 
is governed by a constitutive equation giving the true stress, o(t), when S(t) 
and its past maximum are known. 


Sp(t) = max S (7) 


(52) 
O0<7 «st 


Without loss of generality, this constitutive equation is written in the form: 


s 
0 = Sp (2 : sa) (53) 
n 


and it is assumed that the damage function, F, satisfies the following 
conditions: 


F(1, = 1 
eee (54) 
F(x, Sy)<x for x<1 
These restrictions imply that: 
o(t) < s(t), ie 
also that: 
(56) 


On (t) = Sn (t), 


and, that the following conditions are equivalent: 
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an i) of(t) = S(t) 
Le 4i) S(t) = Sp(t) (57) 
{ iii) o(t) = op (t) 


¥ where o,, is the past stress maximum, defined analogically to Sane The inequality 
" equation (54) asserts that the material actually softens, while equation (56) 
indicates that this softening occurs when and only when S(t) < S,(t) (or 
equivalently o(t)<a,(t)). The results of equation (55) and (56) show that one 
may equally well use the true stress, 0(t), as the damage parameter. 


= Equation (54) and the fact that the first relation of equation (57) implies 
the third are direct consequences of the hypotheses laid down in equation (58). 
To verify equation (55), note that if the maximum of S on the intervalO<;r<t, 

me occurs at 7 =&, then: 


x) Spit) = S(E) (58) 
r) 

ty 

os thus using equation (57) in (52), and recalling that equation (53): 

o(€) = S(E) F(1, Sq) = S(E) (59) 
“ which, by virtue of (57) and the definition of Sp, implies that: 

es a (E) = S(E) = Sq (t)>S (A)DO(A); O<SAKES (60) 


proving equations (54) and (55) and the first two equations of equation (57). 
7 To establish the third relation in equation (57), note that if: 


Sm(t) = S(t) 


which implies that: 


es S,(t) = a(t) 
: 223 
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because of equation (58); then, since 
Sm =%m (61) 
one would have that: 
On (t) = a(t), 
Conversely, if 
o(t) = o,(t), 


then: 


S(t)>0(t) = om(t) = Sa(t) 


so that: 
S(t) = Sat). 


Returning to the constitutive equation (52), it is interesting to consider 
the special case in which the damage function depends only on S/S,: 


sl Sn) . a) (62) 


which is a master damage curve of the type considered in the rate-independent 
model discussed previously. 


When the virgin stress obeys an elastic stress-strain relation: 


(63) 


then: 


224 


eg Og Oa mw ee Leal Slew ee aS im tee ae OR oy, 


a fe ee ec gn ee em EE mR cme Rae oe Ms SS gg eg Ey Fe NT ah ee tae Ne 


Sn = Een 


(64) 
in which €y is the past strain-maximum, and equation (52) yields: 
€ (6 
Os E€n F(e5 ry Een) 5) 
so that, defining: 
€ € 
Fe (= En er(=, Ee] (66) 
n nm 
leads to the starting assumption of Gurtin and Francis (Reference 26): 
Cs re(s, en) E€n (67) 
m 


presented earlier as the rate-independent model. 


Although implicit in equation (52) is the assumption that the functional 
form of F would be the same for unloading conditions as for reloading, it was 
found experimentally that a different damage function is needed for each of 
these processes. Actually, there is more than one way of obtaining the same 
damage function. For TP-H1011, for instance, the following procedure was 
employed. 


Considering the strain history shown in Figure 95, on the loading portion 


we have: t 
: (68 ) 
S(t) = € fon dt for t<T 
° 
hence, S(t) increases monotonically and, by equation (55): 
Sm(t) = S(t) = o(t) (69) 


and, upon unloading, the past maximum of S is the true stress: 
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Om = 0(T), T<t<2T (70) 
| 
Further, by equation (50): 
T t 
sce =12i | f G(t - 7) dr- [ ce-n a (71) 
fe) T 
or, equivalently: 
S(T +t) = G(t) - o(t) (72) 
with: t+T 
ate) =| lf G (1r)dr 
t (73) 
Letting 04(t) and a5(t) denote the 
€ true stress during loading and unloading, 
respectively, with t in o2(t) measured 
a Ss from the time T at which unloading 
t begins, equations (52) and (67) 
T yield the simple formula: 


Figure 95. Strain History Used 
To Characterize the Damage Function 


do(t) [= - 04(t) ais 
s ——— + Om 


on on 
Thus, summarizing, the stress-softening approach to damage is described 


through the following constitutive equation: 


Ss (75) 
o(t) = Sp i ; Sa) 
m 
where t 
de 
S(t) = ic -7)—(1r) adr 
dt 
° 
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and in which the damage function, F, may be determined from saw-tooth tests with 
increasing peak strains and with sufficiently long rest periods between cycles. 
For conditions of reloading, F is given by: 


02 01 (t) 
— =F » om (76) 
on on 


while for wmloading the following form is employed: 


Oo G(t) -%(t) 
—— = Fj ——_———_,, °n, 


on 


(77) 


Typical curves for vunloading and reloading damage functions of TP-H1011 
are included as Figure 96 and 97, respectively. 


Finally, we point out that the use of this approach to predict the 
response of TP-H1011 yielded results that were far more accurate than those 
obtained with any of the other theories in their original form. 


4.2.4.3 Nonlinear Models Based on Maximum Strain. 
A series of constitutive relations based on the past maximum strain have 
been proposed by M. Gurtin. The precursor of these relations took the form: 


t 
de 

o(t) = f G(t - 7) — (1) dr (78) 
dr 


° 


where G represents the relaxation modulus, and the function e was expressed as 
a product of the reduced damage function, F, and the virgin-response function, 
g, in the following way: 


(79) 


with: 
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F(1,€m) = 1 


and 


en = max €(r) 
Og 7T<t 


so that, during virgin response, since: 


€s em 
one had: 
t 
dg(e) dé(r) 
o(t) = G(t - 7) dr 
dé dr 
fo) 
and, by taking: 
k 2K 
g(e) = > Ay ek 
. k= 1 


(80) 


(81) 


(82) 


(83) 


(84) 


the best values for the a;,'s could be determined using least squares and the 


data of all constant-rate tests. 


To characterize the reduced damage function, F, involved the determination 


of a creep function, J, solution of: 


t 
aJ( 7) 
G (t- 7) dr =1 
ar 
Ce) 
and, such that: 
t 
do (rT) 
e(t) = J J (t - 7) dr 
ar 
fe) 
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Thus, taking the reduced damage function, F, in the form: 


F (x, y) = Fy (x) Fo (x,y) 
M- 1 
Fy (x)= xM + S dm (xt - xM) 


m= 1 
P Q-1 


Fo (x,y) = 1+ > bp yP[x - x2 + > Cg (9 - % ] 
q-2 


‘ 


(89) 


ps1 
and equating equations (71) and (75), the coefficients entering F were to be 
determined using least squares and all the saw-tooth data with increasing strain 


peaks. 


When this constitutive law was applied to UTP-19,360B data, it was deemed 
necessary to change the form of the function e, because of the large errors 


observed in the predicted response. 


The last of a sequence of modifications yielded the following stress-strain 


woe fosen eed frgges) fo 


° 


law: 


where, as before, G was the relaxation modulus, and: 


F (1,€m) age | 
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In the present case, the virgin response was given by: 


: . . Fm 
o(t) = { oe ~ 7) K (€m, €m) — 47 (93) 
dr 


° 


while, the damage response, for which €m remains constant, took the form: 
: F 
a(t) = K (em, 0) Git - ne (x,y) 47 an 
fr) 
in which: 
K(e,€) = 1+ Ay (€ -€) + Ap (€ 7 €o)@ + Ag (€ ~ €9)3 + 


: (95) 
+ € (B}€ + Boe€2 + B36 3) 


F(x,y) = a(x) [: + (De y + D7 y?) | - x3 + De (x? - x3)| | (96) 


4 
a(x) = x5 + > Dm (x™ - x5) (97) 
m= 1 
with: 
x = €/E€n (98) 
and: 
y= €m (99) 
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A set of stress predictions obtained for UTP-19,360B, with the resulting 
version of the theory, is included in Figures 98 through 104. Figure 98 
is for the high-to-low dual rate test, Figures 99 through 101 are segments 
of the long-duration similitude test, and Figure 102 through 104 are segments 
of the three step relaxation test. 


Since the dependence of the function K on the strain rate was felt 
to be artificial, the treatment of damage was revised in the manner explained 


next. 


4.2.4.4 Current Model 
The latest version of M. Gurtin's nonlinear stress-strain law is based 


on a strain-dependent relaxation function and has the form: 


t 
o(t) = f sfeco lect - tT) at 
° 


G(E, t) = Gp (t) + Gy (€, t) 


relaxation modulus 


= correction modulus, defined as: 


N 
-t/ 
th= Sn Ce) Ce) : 


n=1 


P 
P 
pet 


For this material, the virgin curve ( €= Em = e) is: 


with o,, the linear viscoelastic stress: 
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Figure 103. 
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o p(t) = J Gr(T) € (t - 7) dt 
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and the correction stress, 0,, given by: 
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Og (t) = J vlecn, r fect - 7) aT 
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Hence, to characterize the Sirdin response, only constant rate tests need 
be employed. In this instance, o and o, are known, so that from equation (88), 
Og may be computed, and equated to equation (90) using the fact that eis a 
constant; i.e.: t 
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which, upon recalling equations (91) and (92) becomes: 
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ey Vn (t) = > Anp€ | e P(A) e da 
p=i1 


° 
bee Furthermore, since for a constant-rate test: 
e(X) = €A (110) 


it follows that: 
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Clearly, equations (89) and (92) to (94) may be used to determine the 
coefficients anp appearing in the definition of the correction modulus. The pro- 
cedure suggested by M. Gurtin to accomplish this is as follows: 
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2. Select the degree, P, of the series expansion of the correction modu- 
lus, as it appears in equation (87). 
3. Use the € test and the approximation: 
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5 Og (t) = Wy (ey, t) 
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to find the ay. 
y's 4, Use the Eo test and the approximation: 


Og (t) =H, (En, t) = W (Eq, t) (116) 


to find the 82p° 
5. Use the €3 test and the approximation: 
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to find the a3p, and so on for the ayp- ++ aANp- 
6. Iterate this procedure if necessary; that is, define: 


N (118) 


so that W is known. For each ramp test, define: 


— 


ef %ery (119) 


and repeat the above procedure using 0, to find constants app * The new 
values of the anp are: 


(Anp) new = Anp + Anp (120) 


7. Repeat the process if necessary. 


It is important to point out here that numerical difficulties may be 
encountered in applying this technique to characterizing the virgin response of 
solid propellants. In fact, some convergence problems were faced in connection 
with the UTP-19,360B data. Moreover, characterization of the damaged response 
calls for a large number of cyclic tests over a wide range of rates. This 
inereases the convergence difficulties. The model was employed with the 
constant-rate tests only for this reason. Figures 105 to 108 show the results 
of the stress predictions obtained with the current version of the model. The 
first two plots correspond, respectively, to the lowest and highest rates 
available in the data base at ambient temperature. 


So far in the program, none of the models developed by Gurtin has taken 
into account the effects of temperature on propellant response. However, the 
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time-temperature superposition principle was tested with the current version of 
the theory. The results appear in Figures 107 and 108 for the thermal tests at 
123 and 40°F, respectively. The use of the superposition principle breaks down 
at a strain of about 39%. This was apparently due to the ambient temperature 
data base being limited to a low strain level. Also, it might be necessary to 
allow the relaxation function (Gq) to depend on the glass transition temperature 
(3) 


Gy = G(€, Tg» t) 


4.2.5 Russian Approach to Physically Nonlinear Viscoelastic Solids 

The Russians have explored two approaches for characterizing damage effects 
in solid propellants (References 7 and 8). They are a general functional 
approach, and a kinetic equation of evolution for damage. Both approaches are 
based on internal-variable concepts, and either approach appears general enough 
to also incorporate cumulative damage and propellant response under multiaxial 
stress states. However, the general functional approach may be of little 
practical engineering value, because a very extensive testing program may be 
required to evaluate material parameters. This approach requires introduction 
of damage measures, which should reflect microstructural damage mechanisms, and 
a damage functional which characterizes the accumulation of damage or defects. 
The damage functional is then expanded into a series of multiple integrals in an 
analogous fashion to that followed by Green and Rivlin for nonlinear materials 
with fading memory. Herein lies the difficulty. Even assuming isotropy, four 
to six different types of multi-axial tests are required to evaluate the 
required material property functions for a first-order theory. Although the 
approach has theoretical merit and may even have some practical application in 
the future, its pursuit was abandoned in favor of the kinetic approach. 


The essential feature of the kinetic approach is to introduce the degree of 
damage into the constitutive equations as a reduced-time parameter in the same 
way that temperature is introduced as a reduced-time parameter for the thermor- 
heologically simple materials in linear thermoviscoelasticity. Damage is then 
defined in terms of some strength parameter of the material, and the degree of 
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damage is characterized through an equation of evolution for damage, as 
explained subsequently. 


4.2.5.1 Original Model 
The one-dimensional constitutive equation taken from the Russian literature 
by W. L. Hufferd as a means of predicting the response of physicaiiv nonlinear 
viscoelastic materials may be expressed by: 
t 


de 
a(t) “| E (t' - 7r') —(1r) dr (121) 


where dt 


o= stress ° 
€ = strain due to mechanically applied stress 
E(t) = relaxation modulus 


E(t) = Eg + Et 


and t 


de 
t' = 7! f 
an(n ()] (122) 


Tt 


represents the damage-reduced time, which is arrived at in the manner described 
next. 


First, a normalized damage function w=w(t) is introduced through the 
following kinetic equation of evolution: 


oe ty eh () P(t) 
dt (123) 


in which it is further assumed that: 


t 
f(t) = J F(t - 7) 6[og(r)]ar (128) 
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together with the conditions: 
w(o) 


w(t) = 


indicating that no damage exists at the initial state, and that failure occurs 
at time t*. 


Next, equation (97) is integrated withw(o) = 0, leading to: 


W t 


dw 
eee fs ve (126) 


h(w) 
° ° 


Setting t = t*, so that w(t®) = 1, and substituting equation (98) for f(7), it 


t# g 
J at [re - 17) 6 [oo (r)lar 
0 fo] 


is obtained that 


1 
dw 
h(w) 
° 
If now the function F(t) is assumed to have a power-law representation: 


Equation (101) can be written in the form: 


S ei hE IO. 
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t* é 
| ak J Fo(— - 7)@¢lo0(7) | dt (129) 
° ce) 


and integrating with respect to &, assuming that the order of integration may be 
interchanged, one arrives at: 


t* 
Fo . 
(t® - r)1 +m¢lo, (7) dt 
1+m (130) 
° 
= 1 
1 
dw 
(% 
° 
which, for the case where 0. and ¢[0. ]=¢@[0. J are constant, becomes: 
Foo (%) (te) 2m 
pny o eae Sere ae (131) 
(1+m) (24m) 1 : 
ao 
h(w) 


where t*, is the time to failure under the constant stress 0,. Thus equation 
(104), in this case, may be written as 


t* 
(132) 
(2 + m) (et cial Pe =1 
ita 
°o 


If the time to failure under a constant stress,o,, has the power-law 
representation: 
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then equation (106) can be put in the form: 


t® 


ei J (t® - 7s o,eta) drs pre (134) 
y 2+m 


° 


So that, motivated by equations (106) and (108), the degree of damage accumulation 
may be introduced through the expression: 


FE: t 
N(t) = (2 + »f (t 


°o 


in which, obviously: 


n (0) 


n (1) 


The function 7(t) can be related to the 


This means that 7 represents the relative damage in the load history for the 
power-law representation of t®,. From equations (107) and (108), equation (109) 
may be written as: 
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<i and finally, the influence of damage is treated as a reduced variable by 


introducing the modified time, t', defined by 


ie dt 

a Ae ee (139) 
4 7 

ort an|nct)| 

wy 

<j on which equation (96) is based, and where the shift function due to damage, 
. an, depends on the material at hand. 

Despite its rather appealing physical and mathematical foundations, the 
| original model did not do any better than Linear Viscoelasticity when it was 
oy used to predict the response of either TP-H1011 or UTP-19,360B. In those 

re instances, a linear expression: 


) anz=1-7 (140) 
: 

Al and an exponential form: 

Be an = e7” 

Bs , 

a were used for the damage shift function. 

Es The partial failure of the Russian approach to reproduce solid-propellant 
“5 behavior made it necessary to change certain aspects of the theory, as explained 
oe next. 

”™ 

¢ Current Model. One revised version of the Russian stress-strain law takes 
ra the form: : 

2 t-r\ de 

‘T4 o(t) = E(—g—]— (1) dr 

~ a, / a (141) 
vg ° 
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i where ah, is a damage-related shift function, assumed to depend only on the 
= current state of strain; specifically: 
E oo ; (12) 
mi azar ett), e(t)| 
a Clearly, if: 
E(t) = E “n (143) 
4! 
oH then equation (113) becomes: 
y 
, t 
eR e on “¢ 144 
al o(t) = Ee e(t) + (a, )” Eo | (¢ - 7) ®—~ (1) ar Cra) 
a dt 
4 ° 
my which resembles the classical approach of the softening function used as 
. a stress-correction factor. 
4 Another revised version of the Russian approach consists of retaining most 
N3 aspects of the original law, but constant strain-rate data are employed to 
4 express the time to failure as: 
* : . 
Ps to (+) (145) 
mh € 
») and equation (108) 1s changed to: 
. (2 + m) (146) 
" 
9 a(t) = | (t= 7) ™O(r) ar 
5 fto 
fd 
a] fe] 
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pS $(t) = e(t) = et 
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Thus, evaluation of equation (117) at t = t® yields: 


tate 


t* 
(2 + m) 


oa J (t® - 1) (rarest 
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which, through a change of variables and after some algebraic manipulations, may 
be integrated to: 


ae 2+m 
(to) =3+m (148) 


The solution for m, as a function of strain rate, is easy to obtain using 
equations (116) and (119) (as presented in Figure 109 for UTP-3001). 


In much the same way, integration of equation (117) for the relative damage 
unetion, n(t), leads to: 


t 3¢m 
n(t) *(—r) (149) 


to 
Hence, using constant strain-rate data to express the time to failure does 
simplify things, but a major assumption would still be needed regarding the form 
of the damage shift functio as it was in the original set of equations. In 
this context, it is important to note that the linear expression: 
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were used for the damage shift function without success. For this reason, the 
modified version used to run the stress predictions included in this report, 
corresponds to equation (115). 


4.2.5.3 Stress Predictions 

Figures 110 to 114 show the comparison between the observed response and 
that calculated using the present theory. As may be seen, the predicted 
response is quite accurate in all cases considered, which include constant- and 
dual-rate tests, as well as a short-duration similitude loading. 


4.2.5.4 Material Characterization 
As may be gathered from equation (115), the simplest version of this theory 
requires the knowledge of only two material-property functions, to wit: 


1. The relaxation modulus, and 
2. The damage shift function: 


a def ° 
(ay) = an = ay [ece), € (t) (152) 


which is determined in the following ways: 


ae From constant strain-rate tests, to correct the stress response during 
loading; 

b. From a relaxation test at a large strain level, to account for 
healing; and 

C. From a constant strain-rate cycle carried out to a large strain level, 
to more adequately reproduce the hysteretic behavior of the 
propellant. 


The damage shift functions corresponding to UTP-19,360B are shown in 
Figures 115 to 119. The first three of these plots represent typical curves of 
a, for low, intermediate and high strain-rate tests, while Figures 118 and 119 
give the correction curves for relaxation and unloading, respectively. 
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Figure 118. Softening Function During Relaxation 
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4.2.6 The Swanson Nonlinear Constitutive Law 
4.2.6.1 Original Model 

The framework for this theory was established by taking into account some 
typical behavior aspects of high-elongation propellants as indicated in Refer- 
ence 9. The principal features considered were: (1) the usual viscoelastic 
dependence of the response on the strain rate, (2) the ability of the solid pro- 
pellant to sustain large strains, (3) the marked deviation of the solid- 
propellant response from that associated with Linear Viscoelasticity, as 
evidenced by the large hysteresis exhibited under cyclic loading by many solid 
propellants, even at small strains, and (4) the dependence of the stress-strain 
response on superimposed pressure. 


Although it is not essential to have done it this way, the capability of 
handling large strains was incorporated into the constitutive equations by using 
the cauchy-stress tensor (9) as a measure of the state of stress at a point. 

Its conjugate, the left Cauchy-Green deformation tensor (B) was used as the 
measure of straining. The Cauchy stresses, defined in terms of force per unit 
deformed area, are also called "true" stresses. In a principal coordinate 
system, B takes om the diagonal form: 


(153) 


in which the A,'s are simply the extension ratios in the principal directions. 


The remaining aspects of the observed response of solid propellant were mod- 
eled through the use of a softening functio as a stress correction factor. The 
major constitutive assumption in this theory relates the second invariants of 
the deviatoric stress and deformation tensors through the equation: 


(154) 


VIIg" = (f) (g) 
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This separable form has been used previously (References 11 and 29) and 
is motivated by the fact that the constant strain rate tensile curves are 
roughly similar. 


In equation (122), f is the following viscoelastic function: 


t 


fe (t ~ 7) MIB dr (155) 


dr 


ba] 
" 


° 


with G being the relaxation modulus in shear, taken in this theory as one third 
the tensile relaxation modulus, and: 


@ = softening function 


O43 - (Gy%/3) 5445 = deviatoric stress tensor 


Q 
= 
Ce 
iy 


Byy - (Byx/3)544 = deviatoric deformation tensor (156) 


o 
~ 
= 

1) 


1/2 
a 2 2 
IIa = -[o 11 A292 + %29 %33 + 233 ayy] + QD + %53 + %34 


Second invariant of tensor a= g, B 


Now, g is a function of deformation and pressure (mean stress) and can be 
considered to be primarily a strain-softening function. It is defined as that 
function of the invariant \ IIg' that will force the viscoelastic Cauchy stress 
to coincide with the experimental results; thus, unloading hysteresis as well as 
the effects of pressure may be readily incorporated into this theory, simply by 
obtaining the corresponding forms of the softening function under such 
conditions. 


The softening function corresponding to virgin loading is obtained by fit- 
ting the model to wiaxial tensile tests at constant cross-head speed. Under 
these conditions, the deviatoric stress invariant reduces to: 
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oe ae (157) 

where, again 0,, is Cauchy stress. 

Assuming incompressibility: 
A, AQAZ = 1 (158) 
and noting that: 

do = A3 (159) 

the deformation invariant becomes 
II,' = —(\y - —) (160) 


Taking the rate of change of this invariant as being approximately 
constant results in: 


t 
r= ¥3x[ —@(t-1) aT et 


° 
so that, from equations (155), (158), and (160), the following is obtained: 
t 


(162) 
ae Vax face ar 
V3 


° 
from which the softening function, g, may be obtained. The assumption leading 


Ka to equation (160), that the time-rate of change of YIIp' is approximately 

be constant, need be guarded against for conditions of changing strain rate. For 
A> example, as in dual-rate tests where viscoelasticity does not predict as fast a 
Me response to the rate change as is experimentally observed. 

ar 
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The modification to linear viscoelasticity necessary to accommodate this 
a behavior is as follows. The response of the function f in equation (163) to a 


‘ constant time rate of change of the deformation invariant is defined as f,. It 
x can be expressed as: 
E é 
: fo =VII vf Gre1 (§ - 7) dT (163) 
i) 


where € = VII pg! /VIIp' 


Atha 


‘ 


The modification to the f function is done in an incremental manner through: 


= fmodified = f +8 [£6 - ANE B' (164) 

s 

bi and the following incremental relationship is used: 

bs dfmod (165) 
‘< + @lmo 

i 7 =f) ——at 

a t + dt t dt 
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we 
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The parameter $ governs the response of the f function under changing strain 
rates. Asf$>0, the response is analogous to linear viscoelasticity. 


oa 
4, 
“a 


The algorithm developed by Herrman and Peterson (Reference 30) has been 
mY) used to implement the calculation of the convolution integral for f. In brief, 
let the shear relaxation modulus be represented by a Prony series as 


sp s 


“ m 
* (166) 
i G(t)= > gerast 
i=1 
a4 Then let the f function at time t, be given by: 
Z tn 
3 S @a5 (ty - 7) ayii B' (167) 
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A recursion relation can be easily developed to compute f(t,) (Reference 30). 


Let a 
iis 1 
and 
tn 
Tad | Gre" az(ty - 7) 8 Vir Bt a (169) 
° a7 
then 


-O, At : Gy - a sAt 
Inaee + 7 Intai+ Vary af i. © " 
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giving for the change in these terms 


. Gy -a@ =; ¢ 
4Iy 41 = VII Bh e [ -e a + In-1,4 + . sen = ] (171) 


a 


* 
- 
a Fa 


mn which is directly analogous to linear viscoelasticity. The modification 
oe proposed above can then be implemented as 

bs (Aly 1) modified = 41n,i +6 [Ton,4 7 Tp,a| A¥FTp' (172) 
ao, 
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is and the I terms can be incremented according to: 
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oy (173) 
o In,i = In-1,i + (4%n,i) modified 

a This has the effect of changing each term of the series so that it approaches 
< the value it would have been if it was always at the new strain rate. Note 
ee again that (as discussed in Reference 30) varying temperatures can be 


incorporated into the time scale as usual. 


Unloading tests required a further refinement of the model. For lack of 
more detailed information, the parameter 8 may be taken as zero for wnloading 
Me states (i.e., states in which VIIp* is decreasing). The large amount of 
- . hysteresis seen in load-unload cycles is then modeled in part by the hysteresis 
inherent in linear viscoelasticity primarily through the g function. This is 
accommplished by giving g a different value when the deformation invariant YIIp' 
SS is less than its maximum previously achieved during the loading history. If 
V15"nax is the current maximum value, the function: 


‘ Pat on 
Poe or 
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-Re 
< 
* 


BS ge = @(VII p' max) ) - Cy [; -VII B' /NIlp: = || 
aoe (174) 
eis 
, provides plasticity-like behavior. 
on 
; Py The behavior of the g function for unloading and reloading conditions is 
Se illustrated in Figure 120 (taken from Reference 9). 
The Swanson approach (Reference 9) 
i was used with an only limited degree 
“ S of success to predict the stress 
va, c 
ey 2 response of TP-H1011 and UTP-19, 360B 
i= 
= under several strain histories. 

a F 

a = 

oa 8 In the case of TP--H1011, the errors 
wT £ 

“i § in the predictions were believed 

4 
. to be due to uncertainties in the 
a 9 value of the changing-rate coeffi- 
os 0 0.4 08 12 16 2.0 lent (8) Th dat il 
Me Second invariant of B’ (v1 B) oe . Chena ee eae ee 
wa Figure 120. Effect of Deformation able to determine @ directly for 
= and Pressure on the Strain this propellant. 
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oe It was possible to characterize UTP-19,360B in a complete fashion. The 
corresponding predictions were not any better than those obtained for TP-H1011. 
This led to changing the law as discussed below. 


~ 
s 4.2.6.2 Current Model 
in Analysis of the stress predictions, carried out for UTP-19,360B with 

sa the original Swanson theory, revealed the importance of several inadequacies and 
x oversimplifications listed below. 
1. The softening function (g) should depend not only on the strain and 
uN pressure but also on the strain rate. 
i 2. The softening function, as defined by equations (155) and (163), 

~ should be different for unloading than for reloading. 

am 3. The softening function for unloading or for reloading should never 
Ly become zero for conditions of tensile straining only. A zero value 
. could occur with the softening function defined by equation (174). 
3 4, The healing process observed during relaxation in solid propellants 
a like UTP-19,360B was not taken into consideration by the original 

“ Swanson theory. 

S 5. The reverse-recovery observed in solid propellants during relaxation 
“5 or rest periods that follow an unloading process, only poorly modeled 
4 by classical viscoelasticity, is not considered in the approach 

ne by Swanson. 

“2 6. The changing-rate coefficient (8) is more a mathematical device than 
xs it is a material property. If the softening function is made to 

depend on the strain rate then 6 need not be used. 

. 7. The use of a softening function as a stress correction factor 

Sy eliminates the need of using the Cauchy stress (0) and the nonlinear 
o measure of stretching (B). 
a All these observations were incorporated into the original stress-strain 
aa law but the general form of the corresponding equations remained the same, 

ji namely; 
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° 


valid for one-dimensional loading, with: 


t (176) 

at | 
St - Sr = Tee et 
2 ar|(t)| 


representing temperature-reduced time; and where the time-temperature shift 
function was taken in the power-law form: 


(= =) (177) 
ary =| ——— 
T-T, 


in which Tp is the shift reference temperature, T, and m are material 
parameters, and T is the current temperature. 


The modified version of the Swanson theory was most successfully used 
to predict the response of UTP-19,360B, as explained next. 


4.2.6.3 Stress Predictions 

The degree of accuracy of the predictions made with the current version 
of the Swanson approach may be realized by examining Figures 121 through 
130. The first two figures correspond to the lowest and highest constant- 
rate tests available. Figure 123 shows the predictions corresponding to 
a saw-tooth test at constant rate and increasing peak strains. Figures 124 
and 125 present the results for the dual-rate tests while Figures 126 and 
128 pertain to the long- and short-duration similitude tests, respectively. 
Figure 127 shows the three-step relaxation test. Finally, Figures 129 and 
130 show the results obtained for constant rate tests at 123°F and 40°. 

(Text continued on page 285.) 
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RU & - 


Pee 


Ss 
oy 
ad 
vy They testify to the fact that the time-temperature superimposition principal 
+ may be used without sacrificing more accuracy than is already lost in fitting 
e. equation (143) to the very limited time-temperature shift data. 
Pa 
Pa) 
cs 4.2.6.4 Material Characterization 
es According to this theory, only the following listed properties are needed 
ae to characterize a solid propellant completely: 
fe 
oe 1. The relaxation function, G, defined as: 
, Ere1 (t) 
P” Gc (t) = —— (178) 
a 3 
Re where Enel (t) is the linear viscoelastic relaxation modulus. 
2. The softening function, g, defined and obtained as follows: 
, a. For loading conditions:. 
ral 
e g = 8 (e, 2) (179) 
8 L 

and it is obtained from a sequence of constant-rate tests at at least 
= three different rates that span the range expected in the applications. 
. b. For unloading conditions: 
eo Bu = & (€/€max) (180) 
3 
sti 
ca 

where *,,, represents the maximum strain previously achieved 

io during the loading history. The g is determined from the 
3 unloading portion of a loading-unloading cycle carried up to an 
4 intermediate strain level. 
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Ce For relaxation conditions: 


Br = Br (t - to) 
(181) 


in which t. is the time at which the relaxation process begins 
amd g is evaluated from a relaxation test at an intermediate 
strain level. 


In addition, for relaxation after partial unloading or during rest periods 
starting at t = t,: 


mn 1 (182) 
Gr = 
Gr (t = to) 


Also, the stress-correction function for reloading is taken as a linear 
function of the relative strain. It is a straight line from the point where 
reloading starts to the point of maximum loading over the past history. 


4.2.6.5 Addendum to Swanson Theory 


Three-Dimensional Version of the Model 


The (general) constitutive assumption used to relate the deviatoric compo- 
nents of the stress and deformation tensors, takes the following form: 


O45 Bs, 
Vit, *Vitp i, j=, 2, 3 (183) 
together with: 


Vitj= (g) (f) 
(184) 
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- 7. a ee 
el a 


or, equivalently: 


t 


Vit = (g) { G(t -7) aVIIpr at (185) 
ar 


° 


O45 = i-j component of the deviatoric Cauchy stress tensor 
BY; = i-j component of the deviatoric Left Cauchy-Green deformation tensor 


x ,, IIp, = second invariants of the deviatoric stress and deformation tensors 
x 

2 with: . 

an ' 1 

ven O45 = %%4 +— (944 + %99 + 933) 55 4; i, j= 1, 2, 3 

- ? (186) 
= 


| def ’ . 
IIe = ~[o11 22 #5 933 + O22 033 | + (842)2 + (843)2 + (053)? 187) 


and similarly for Byj and I1g'; and in which: 


ne def}1 for i= j 
Si 545 = 
ne O for i¥j 


also: 


Ss & = softening function that depends primarily on the strain level, the strain 


rate, and the applied pressure. 


and: 


; ect) o* e¢ty/3 
(189) 


where E(t) represents the tensile relaxation modulus at a small strain. 


Og eee Dealt ast wi ain? ak. Yee ea ie (Set ee tomey om Oil Bie lat sete Sagi LG RE a Me ait ae 
Rat She re a 2 et ee Bt Sg al eS a aS sg ete ek ie Sh SR Sg ed 
a ee eiea ie Boe? BP Fa WO Se cre a ee gn Nar ek Nee a ar ne a ne ae a ie 
Fe te Mi ca Fan wy Ma ee Ty Te ME Se Bo ce OXe Se eee 


Lar a ta ea Ga ae ae oe en ee, Pa 


ae According to the constitutive assumptions (183) and (185), the distortional 
| behavior of the material is completely characterized through the softening 

ie function, g, and the relaxation function, G, which may be evaluated from one- 

n dimensional tests, as explained in the previous section. Indeed, the stress- 
strain relations set forth in equation, (183) and (185), reduce, as they should, 
to those employed in the one-dimensional version of the model. 


To complete the theory, an assumption is still needed about volumetric 


: behavior; and although time-dependent bulk response may be important in some 
rf applications, the elastic relation: 


1 1 
oon +959 + 033) = K ([AyA2A3| - 1) (190) i 


may be employed; in which K is the bulk modulus and the A,'s are the stretch 


| 
os ratios. | 
| 


~ For an incompressible material (and solid propellants are nearly 
inec apressible): 


Cs, (191) 
A Aq, A2A3 21 


; so that equation (190) breaks down, and the stress tensor has to be considered a 
on function of the mean pressure (0,; + [%po + 933)/3, as well as of the deforma- 
tion tensor, leading eventually, to a stress-strain law of the form given in 

o equation (183). 


ms Application of the Model to Two-Dimensional Problems 
45 In order to use the stress-strain law presented in the foregoing section, 
as one must have available the deformation tensor at each point of the continuum 


where the stresses are desired. This solution in terms of deformation may be 


arrived at numerically through, say, finite elements, or, analytically. 


o. The accuracy with which the present constitutive theory may predict the 
e two-dimensional response of solid propellants may be seen in Figures 131 to 
a 

s 
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136, which correspond to constant strain-rate tests of strip-biaxial samples of 
UTP-19,360B. The first three figures belong to tests performed at a nominal 
crosshead displacement rate of 0.02 in./min. at 40°F, 70°F, and 120°F, respec- 
tively; while Figures 134 to 136 show the results for a crosshead displace- 

ment rate of 0.2 in./min. at the same low, intermediate and high temperatures of 
40°, 70°, and 120°F. The plotted data refer to the direction of applied loading, 
which is also the direction of maximum principal stress (and strain). The 
geometry of the strip-biaxial sample used is as presented in the following sketch. 


(direction of loading) 


050/011 = o> 0.481 (192) 


0 MD ME OE A a a A a A 


2b 


ae, 


€20/€ 41 =O, = 0.025 


The stress- and strain-axiality factors, ¢g and ¢_, were taken from 
Reference (31), and are valid at the center of the sample for small strains 


only. 


The constitutive relations given in (183) and (185) yield: 


with summation implied over repeated indices. 


Now, under conditions of plane stress of an incompressible material, 
and along the principal directions, one has: 


0 
(20 2 - 04) 0 
0 - (0, +05) 


0 
2 2: <2 
(242 - AZ = Aq 290 2 2 
(2A3 = ry - 5) (196) 


t ° , 7 e ¢ 
II B’ = 7By1 Boo - By B33 - Boo B33 


To evaluate (194), we first write it in umabridged notation, noting that in 
this case, if i # Jj then By, = 0; thus: 
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and using (197): 
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where, from (196) one has, for instance, that: 


20 2 2(aa, hae a, 23 ) (200) 
at 3 dt dt dt 
and similarly for the derivatives of Boo and B33. 
In the previous derivations, the stretch ratios are computed as: 
Ay = 14 €, (t) 
Ape ten (t) 2146, €4 (t) 
Az = I/(AqA2) oy) 


in which ¢€,(t) is the strain 


history imposed on the sample along coordinate 


1, and the last expression of (201) follows from the incompressibility condition, 


equation (191). 


Hence, the first component-equation of the constitutive relation (193) 


yields: 


, 1 
"Ue tg (20 4 @2):8 


or,in view of (192): 
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Using equations (192), (196), and (199) to (202), we obtained the response 
of the biaxial sample in the direction of the applied loading. The plots 


ineluded in this report, show the engineering stress: 
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rather than the Cauchy stress, o. 


4.2.7 M. Quinlan's Theory of Materials with Variable Bonding 


4.2.7.1 Original Model 


ox In developing a mathematical framework for his stress-strain law, Quinlan, 
; in Reference 4, reasoned that since propellants consist of minute rigid particles 


f embedded in a polymer matrix, such materials would respond to a deformation proc- 
_ess with a change in the amount of species to species bonding. He thus proposed 
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to correct the deficiencies of fading-memory type theories by introducing a 
correction term that accounted for the changes in the state of bonding that are 


induced by a def ormation process. His constitutive model then took the form: 


pati 


a” 


ee ae, 


rit 


(204 ) 
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Oo =0¢ + Op 


x in which: 

ES 0 = current stress 

S Of = fading-memory type stress 

e Op = stress correction due to change in the state of bonding 


Motivated to some extent by reaction-rate theory, Quinlan modeled the 
evolution of the bonding state through the following ordinary differential 
equation: 


Pit ig 


. 
aa 


st 


tts 
. 


hea d-u[1 - ee (@ - »)] 


subject to the initial condition: 
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in which 7 represents the state of bonding; a, yw and pv are material parameters, 
4 and 


(207) 
1+e 


is the stretch ratio; with e, the strain. 


The unique solution of (204) may be readily obtained for piecewise linear 


stretch histories, as reported in Reference 4, 


Taking a linear viscoelastic relation for d¢, and considering the stress 
as correction term, op, as proportional to the state of bonding, Quinlan arrived at 
the following stress-strain law: 


t 


“ay o(t) = J G(t - 7) $ (T) dT + Bo m(t) (208) 
~ with: 2 


G(t) = Eo t7n 


n (t) =a |>-» [3 ~er(o- 7) | 


eo 


and 
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m(t 2-0) = 1 (211) 


wy The six parameters: Eo, n, Bo, a, handy, needed in this theory to 

rs, characterize a propellant, may be obtained by fitting the model to the observed 
y response of the material when subjected to a saw-tooth strain history that has 
= increasing peak strains and sufficiently long rest periods between cycles. 

i Alternatively, the studies reported in the literature on the effects of employing 
3 different data bases for characterization, show that the test history should pri- 
8 marily include the maximum expected strain level, the expected range of strain 


rates, as well as rest and relaxation periods. 
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The present model was used to predict the response of TP-H1011 under sev- 
eral loading histories; and it reproduced, somewhat accurately, the general 


trend of solid propellant behavior. 


In an attempt to include healing effects, the underlying assumption for the 
evolution of damage were revised, as explained next. 


4.2.7.2 Current Model 
The theory developed by Quinlan has undergone several changes; mainly in 


the expression defining the evolution of damage. Thus, the form: 
t 


o(t) = | H (t - 1) o(r) dr +C nt) (212) 


) 
has been retained in all versions of this theory, but the rate mechanism 


underlying damage: 
m= P(n,6,0) 1223) 
for which: 


m(d21) 21 (214) 


was assumed to contain a neutral rate, ¢, at which damage remains constant, 


i.e., 


P(7,¢,5) 20 (215) 


This concept then allowed introducing the notion that at rates higher than ¢ , 
bond breakage would take place, while at rates smaller than the neutral rate, 


bond formation would ensue. The details of the derivations leading to the 
specific form of equation (213) are presented next. 
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Equation (215) for the neutral rate, {, may be rewritten as: 
C= Q (m, ¢) 
(216) 
$ = 0 if and only if 7=9 
which, upon expansion in Taylor series, becomes: 


$= Q (1, 7 -6) +29 1-0) (mn - 6) + Q (|7 - 6/2) 
$=0 a(t -¢) g= 0 


The first term on the right-hand side of this equation vanishes by virtue 
of (216), so that, neglecting the higher order terms, and defining: 


def 0Q (7 anes 
usc —_ + + 
O(7 -¢) @=0 


(217) ! 
leads to the following first-order expression for the neutral rate: | 
56 =Hu(n - $$) (218) 
where for bond breakage: 
¢ -~S>0 (219) 
while for bond formation: 
$ - $<0 (220) 
In addition, equation (213) may be cast in the following form: 
(221) 


m= P (1,¢, ¢) = R(s, u) 
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with § given by (218). | 
he,*. 

. Now, under the assumption that the process of bond formation is slower than 
a that of bond breakage, the function R defining the evolution of damage must be | 
a of the form shown in the sketch below, in which the parameter a would represent 
the maximum rate of bond formation. 
no 
ria 
ve 

” . : 
a =@o-$ 
oN 
on Hence, R may be defined through the following differential equation: 

s 
3) 

iy dR(s, u) [Rc ) 

ae =v s - a 

| du saa) (224) 
a 
4%. 
ee where a and vy are positive constants. 

My 

i Integration of (224) yields: 

4 


R(s, u) = a(e’U - 1) (225) 
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Finally, putting (218), (223), and (225) into (221), results in: 
ts afem[ys suv @ -m]-1} (226) 


Equations (212) and (226) subject to (214) were used by Quinlan in several 
ways to characterize the response of UTP-19,360B. One such stress-strain law 
took the following form: 


t 


. (227) 
a(t) = Bo e(t) + [By + Ep co) f (b- nm er) arece | 


° 


in which Eo, E,, Eo, n and C are constants and ¢ is the strain. Although some 
aspects of propellant behavior were better modeled than with the original ver- 


sion of the theory, others were not, and further revisions were necessary. In 
the latest version of his constitutive law, Quinlan used a Prony series to repre- 


sent the relaxation function and changed strain for stretch in the original equa- 
tion of evolution for damage, so that, in summary, the current model looks as 
follows: 


O(t) = oy(t) + % p(t) 
t 


Oy(t) = Ge E(t) + J Gt -7) €(7) dr 


° 
n 
G(t) = S abe 
i=] 
0 p(t) Bev(é -™) €(t) 
erea fore) 4 ;™(0) = 0 (228) 


This final version of the theory has been employed to characterize UTP-19, 360B, 
but has not been used by Quinlan to predict the response of the propellant under 
any loading history other than the characterization test. 
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4.3 CONCLUSIONS 
The best nonlinear constitutive theories available for modeling solid pro- 


. 
ae 
ve 


pellant response were selected. Each of these theories was able to model some 
simple constant strain rate test behavior during the early part of the progran, 
. but generally gave poor correlation with the long time complex load histories 
characteristic of rocket motors. 


A broad nonlinear data base was developed with two solid propellants. This | 
data was used both at CSD and at some of the University subcontractors facili- 
ties. This data base can be used for evaluating any future nonlinear constitu- 
tive theories. All data was collected on an HP-9825 computer and is available 
in a format for a VAX computer systen. 
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This data was used to evaluate and further develop the nonlinear theories 
< for solid propellant response. Each of the theories has been extensively 
modified and now fits both simple and complex uniaxial load histories. 


These theories will be further developed for multiaxial load histories in 
4 the last phase of the program. 
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Appendix A 
MULTISTATION AUTOMATED DATA REDUCTION 


INTRODUCTION 

Automated handling of multistation tester data is accomplished with a 
system of interactive programs on the HP 9825 desk top computer (Figure A-137). 
These programs include data acquisition, stress relaxation-master modulus, 
straining while cooling or heating, straining to failure, and complex histories. 
The acquisition of data and test control are functions of the data acquisition 
program which supplies data to the data reduction programs. The reduction 
programs reduce and output data for a particular type of test history. In 
addition, terminal emulating software for the HP 9825 provide a data link for 
the transfer of data to the VAX-11 mainframe computer. This makes the data 
directly available to the nonlinear constitutive theory programs. 


SYSTEM INSTRUMENT CONFIGURATION 

The multistation data acquisition instruments are configured to provide 
load, crosshead position, temperature and elapsed time data to the data 
acquisition program. The system consists of a Hewlett Packard 9825 desk top 
computer, 3455 digital voltmeter, 3495 scanner, 98035 programmable clock and 
9885 flexable disk drive (Figure A-138). 


The HP 9825 and data acquisition program act as the system controller. The 
controller processes incoming test data and crosshead information and responds 
by sending instructions to other instruments in the system over an HP-IB inter- 
face. Output signals from the tester's load transducers, linear potentiometer, 
and analog thermometers are input into the scanner's programmable relay cards. 
The scanner's relays under command of the controller can be opened independently 
to route output data signals individually to the digital voltmeter where they 
are digitized and read by the program. Crosshead control information from 
output lines connected to the tester's motor-clutch assembly, is supplied to the 
program through the scanner in the same way as the data output signals. These 
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Figure A-137. 


28842 
signals enable the controller to react to changes in crosshead movement and 
direction without relying on operator intervention. The programmable clock 
connects directly to one of the computer's I/0 ports. It provides the program 
with elapsed time data and a program interupt capability for controlling the 
rate at which data is taken. The flexible disk drive provides a mass storage 
medium where data is stored during testing for latter access by one of the data 
reduction programs. 
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a pee PROGRAMS 

es Data Acquisition: 

: z As previously stated, the data 

i Transducers acquisition program is used to collect 

ri data from the mltistation tester. 

et aes Operation of the program involves steps 

a thermometer to initialize the program, calibrate 

F vottmeter the system, and collect and store test 

a Vineak data. 

rans potentiometer 

- Initialization of the program is 

Ne Clutch-motor accomplished by operator entered 

4 aESEOY Figure A-138. information used to identify the 

Sa 28843 ~=particular test and define samples 

being tested. In response to prompts from the computer display, the operator 

re inputs descriptions on test material, crosshead rates, strain levels, and 

ce temperature levels of the test history. Data input on the test samples include 

~ their number, gage length, and individual cross sectional areas, along with 
their channel locations. In addition, the operator enters pairs of crosshead 

os rates and delta strains for each test interval used to compute sampling rates. 

in The cperator also determines how data is taken during relaxation cycles by 

ay specifying whether sampling is to be done in a fixed or log time interval. 

in 

ae Calibration of the system is done by an operator-interactive procedure to 

o ' determine the tester's transducers and potentiometer sensitivites. This 

om involves the operator queuing the program to take readings from the transducers 

ns at differing load conditions. By comparing the change in output signals for a 

ms known change in load, the lb/volt sensitivity of each transducer may be 

2 determined. Similarly, by moving the linear potentiometer probe a known 

distance its in./volt sensitivity is determined. The analog thermometers are 

= not calibrated at the time of testing. These wits output a 10 mv/°F repre- 

oy sentation of the test chamber and internal sample temperature. Calibration on 

3 them is done periodically by the CSD electronics laboratory. For short time and 

isothermal tests, calibration is done once before testing begins. For tests 
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a 

ax lasting over a long time period, a second calibration is done when testing is 
rs complete to enable compensation for drift in the tester's electronics. Since 

: the load transducers are temperature sensitive, for thermal tests two cali- 

~ brations must be performed at differing temperatures to determine the change 

mM in sensitivity per degree change in temperature. 

au 

: When calibration is complete, the program stops operation until testing is 
: ready to begin. On a cue from the operator, zero load and position data is 
a taken and the system's instruments programmed to their initial conditions. The 
2 clock interupt period is set for a sampling rate determined from the initial 

. crosshead rate and delta strain information. Scanner relays are also arranged 
we to monitor the tester's break input voltage. 

: 

3 The program monitors the brake voltage until detecting the brake has 
disengaged which signifies crosshead motion. The clock's counter and interupt 
4 units are then started. Interupt signals are output by the clock at the set 

% sampling rate until changed by the program at the end of the straining interval. 

When interupt instructions are received from the clock, program operation 
¥ branches to a data collection subroutine. The voltmeter and scanner are set to 
MS read output signals from each of the transducer, potentiometer and thermometer 
8 channels. Fifty milliseconds are required to read each channel. Elasped time, 
. read from the clock counter, is taken as the mean time over which the data set 
; was read. The test data is retained in a memory buffer until transferred to a 
Si disk storage file. 

The program continues to monitor the crosshead break and clock information 
channels throughout the test. When a change in crosshead motion is detected, 
2 the clock interupt is stopped. From the test description data corresponding to 
¥ the test interval, a sampling rate is determined and the cle” eset. For log 
a time interval samplings, the clock interupt unit is stopped aiver each reading 
a and the time doubled. 
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Up to 600 data sets may be retained in the computer's memory at me time. 
a Data is transferred to the disk either between data set samplings, if time 
a allows, or when testing is complete. 
= Data Reduction Programs: 
<2 The reduction programs reduce and output multistation data pertinent to 
w particular types of test histories. The programs are stress relaxation-master 
4 modulus, straining while cooling or heating, straining to failure and complex 
Re histories. A description of the strain and temperature histories relevant to 
3 each is listed in Table A-35. 
as 
a Test identification, calibration, load, sample extension, thermal, and time 
a data are supplied to the programs from the acquisition data files or entered 
me directly by the operator. In addition, relaxation cathetometer strain 
ra measurements and the thermal expansion coefficient may be optionally entered. 
= 
a Each program reduces stress, strain, modulus, temperature and elapsed time 
x data when applicable. The method by which each is determined depends on the 
3 test history and amount of information available to the program. 
Rs Calibration sensitivity (S) of the transducer and potentiometer are 
oe determined in general by 
a S = loed/(load output-zero load output) 
ao 
. 
; where load is the transducer calibration weight or potentiometer probe 
displacement. For tests where multiple calibrations were performed, the 
: sensitivities at time t are corrected for electrical drift and thermal 
. 
rs variations with the linear relationships 
ph; 
a (Spinai ~ Sinitiar? 
= Xfinal ~ “initial 
~ 
Bi 
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ES TABLE A- 35 


( Program Test Strain History Output 


‘ 1. Stress relaxation Tabular - time, modulus 


Graphic - modulus vs time 


2. Straining while 


Tabular - time, strain, 
cooling or heating 


temperature, 
stress 

Graphic - stress vs time 
and temperature 


Test History - 
Tabular - time, strain and 
stress 


Graphic - stress vs time and 
straining 


Mech Properties 


Tabular - initial modulus, 
maximum stress 


and strain - 


corrected stress 
and strain - 


rupture strain 


4. Complex histories Combination of strain- Tabular - time, strain 
ing, relaxation and temperature and 
e temperature intervals stress 
an Graphic - stress vs time 
5 and strain and 
2 temperature vs 
time 
where X is time for isothermal and temperature for nonisothermal tests. 
oe Temperature corrections are not made on the potentiometer sensitivity since it 
a is located outside the enviromental test chamber. 
5 Zero load outputs (ZO) for the transducers are also corrected for 
4 
m electrical drift and thermal variation by the same method as the sensitivities. 
t=, 
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Corrections on the zero position output of the potentiometer cannot be made 
since the crosshead can't be accurately returned to its initial position. 


Sample stresses (0) at time t are calculated by 

o(t) = [ transducer output - 20 (t)] S(t)/cross-sectional area 
Sample strain (¢) at time t is determined by 

e(t) = [pot output - 20 | S(t)/gage length 


For nonisothermal tests .a strain correction may be applied using the thermal 


expansion coefficient (a). In this case, the total sample strain becomes 

e€(t) = mechanical strain + [mt) - Tinitial | a/gage length 
An additional correction for effective gage length may be made using cathe- 
tometer measurements of actual sample strains. The correction factor is 
determined as the ratio of the mean intervals in the test history. For 
histories where multiple cathetometer measurements were made, the correction 
factors are linearized to measured strain between them. 

Relaxation and secant modulus (E) at time t is determined by 

B(t) = a(t) [1+ e(t) se(t)] 


where e€(t) is held constant over relaxation test intervals. 


Temperature is ~educed from analog thermometer readings by converting the 
millivolt output to volts. 


Elapsed time is calculated as the difference between when the data was 
taken and when initial loading occurred (t,.) since loading times may vary from 
sample to sample, t,. is approximated by the time of initial straining. 


Once data is reduced, a tabular and graphic summary of the test is output 
by each program. A description of the outputs is listed in Table A-35. To 
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s retain data for future reference and reuse, identification, calibration and raw test 
e data are stored on permanent diskette data files. 

ne Terminal Emulator 

in The program is used to transfer data between the HP 9825 and VAX-11, desk 

~j top and mainframe computers. 

a A link is created between the computer types utilizing the VAX-11's dial- 


= in lines and an RS-232 interface which connects the 9825 to an acoustic coupler. 
The emulator software then supplies the capability of using the 9825 as an 
intelligent terminal through which data may be read from the flexible disks and 
4 sent over phone lines to the VAX. 


Data transfer is accomplished with a VAX program which reads data sent from 
the terminal and retains it in storage files for access by the nonlinear theory 
ee programs. 
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SYMBOLS 


microcrack growth rate shift factor. 
temperature dependent material function 


expansion coefficients of bulk stress in terms of 
octahedral strains 


constant 

initial area 

temperature shift factor (a7) 

(AHETA) damage related shift function 

half sample width 

softening function 

constant 

expansion coefficients of correction modulus 
constants 


biaxial sample width 


Cauchy-green deformation tensor 
bulk modulus and a constant 
deviatoric deformation tensor 


biaxial sample height (gage length) 


softening function 
Chemical Systems Division 
rehealing parameter 
constant 


constants 
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constant 


constants 


modulus 
equilibrium modulus 


reference modulus and normalized coefficient for 
modulus 


relaxation modulus 

E(t), Eng) (t) linear viscoelastic relaxation modulus 
product of F and virgin response function g 
deviatoric strain tensor 


linear viscoelastic modulus 


force 


F(t) . constant rate modulus 


F(E, €m) damage curve at €, damage level 


material parameter 

deformation function 

constant time rate of change of deformation invariant 
viscoelastic type function 

degrees Fahrenheit 

damage function or softening function 


strain magnification factor 


shear modulus 
gage length 
corrected sodulus 


relaxation sodulus 
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shear relaxation modulus 
shear relaxation modulus 
virgin response function 
function of 


strain softening function 


function of damage in kinetic equation of evolution 


Hydroxy Terminated Polybutadiene 


volume dilatation 
contribution to stress at time tn 
octahedral shear strain 


Lp norm 


ereep function 


Joint Army Navy NASA Air Force 


constant 
stress intensity factor 


rehealing parameter 


constant 
constant 


constant 
constants 
constant 
millivolts 


material parameters 
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; N number of cycles 
, n constant 
4 
: P terms of equation under summation 
: P hydrostatic pressure 
1 PBAN Polybutadiene Acrylonitrile 
3 P constant 
Po, Py constants 
: Pas used to normalize Y3 to 1 
Q terms of equation under summation 
: q constant 
RH relative humidity 
8 virgin stress and damage parameter constant 
SQ damage parameter 
: So constant 
Ss. certain measure of damage 
S_-Sr temperature shifted time 
Sy constant 
. JS” damage parameter 
: T temperature 
T peak stress time 
Te material property and shift temperature 
‘ To temperature at t = 0 
: Tp reference temperature 
; 319 
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¢ | 
7 
| 
t 1 
i 
N t time 
i 
ci ty time 
ts time to failure under constant load 
: UTP United Technologies Propellant 
x = €/Een 
x width position from center of biaxial sheet 
X, root of ¥3 
y BCiy 
: Yq. Yq, Yo, ¥3 functions related to damage 
a coefficient of thermal expansion 
a constant 
B constant 
B parameter 
Y shear strain 
AL change in length 
{ 544 Kronecker delta 
€ strain 
if 
; €4 principal strain 
: €> lateral strain 
i € 414» €22, €33 principal strains 
: €o pseudo strain 
; 
‘ € ns “max maximum strain 
2 € (t) strain at time t 
€y strain of unfilled polymer 
t 
; €g strain due to mechanical stress 
; 
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n compressibility (7 = 0 is incompressibility) 
n(t) related to damage function, 

r extension ratio (1 +e) 

healing correction factor 


r softening function 


é 
‘d 
J 
4 
P 
‘ 
4 A 


r width to height ratio 
, g reduced time 
g ratio of position to half sample 
3 width (x/a) 
Lo second invariant of tension (a = a, B) 
S summation 
Co] engineering stress or stress 
0 Cauchy-stress tensor 
OB bulk stress 
G4 stress correction 
og j deviatoric stress 
‘o' deviatoric stress 
A) linear viscoelastic stress | 
06 constant stress 
‘ Ck bulk stress 
4 Or linear viscoelastic stress 
| o(t) stress at time t 
0 ¢(t) fading memory stress 
t reduced time 
T shear stress 
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shifted tine 

function of loading 

nth component of stress correcticn 
normalized damage function 

second invariant cf tensor (a = 0, B) 
second invariant of deformation tensor 
seoond invariant of deviatoric stress 


denotes absolute value 
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